Draft version February 1, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 03/07/07 



THE TOOMRE SEQUENCE REVISITED WITH iTST/NICMOS: NUCLEAR BRIGHTNESS PROFILES AND 

COLORS OF INTERACTING AND MERGING GALAXIES 1 

Jorn Rossa 2 ' 3 , Seppo Laine 4 , Roeland P. van der Marel 2 , J. Christopher Mihos 5 , 
John E. Hibbard 6 , Torsten Boker 7 , and Ann I. Zabludoff 8 

Draft version February 1, 2008 

ABSTRACT 

We discuss the near-infrared properties of the nuclei in the 11 merging galaxies of the Toomre 
sequence, based on high spatial resolution J, H, and K imaging data using NICMOS onboard the 
Hubble Space Telescope (HST). The observations are less affected by dust extinction than our previous 
HST/WFPC2 observations and offer higher spatial resolution than existing ground-based near-IR 
data. Nuclear positions are generally found to be consistent with those reported from data in other 
wavebands. In NGC 7764A, we detect for the first time two nuclei with a separation of about 260 pc, 
consistent with its placement in the middle of the merging sequence. We see a marginal trend for the 
nuclei to become bluer with advancing merger stage, which we attribute to a dispersal of dust at late 
times in the merging process. Our data also indicate a statistically significant trend for the nuclei in 
the sequence to become more luminous, within an aperture of fixed physical size and after correcting 
for dust extinction, with advancing merger stage. We derive -fT-band surface brightness profiles for 
those nuclei for which the morphology allows a meaningful isophotal analysis, and fit the profiles with 
a "Nuker law" for comparison with other samples of galaxies observed with HST. The majority of 
the nuclei have steep profiles that can be characterized as power-law type. In general, the Toomre 
sequence galaxies tend to have steeper profiles and higher central luminosity surface densities than 
E/SO's. Our findings can be explained if the Toomre sequence galaxies have newly formed stars that 
are concentrated toward their centers. We derive V — K color profiles for the nuclei to further address 
this possibility, but find that the large amounts of dust extinction complicate their interpretation. 
Overall, our results are consistent with the generic predictions of A-body simulations of spiral galaxy 
mergers. If left to evolve and fade for several Gyrs, it is possible that the properties of the Toomre 
sequence nuclei would resemble the properties of the nuclei of normal E/S0 galaxies. Our results 
therefore support the view that mergers of spiral galaxies can lead to the formation of early-type 
galaxies. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: interactions — galaxies: 
nuclei — galaxies: spiral 



1. INTRODUCTION 

It is commonly assumed that the merging of two 
spiral galaxies following a close interaction will even- 
tually result in their transformation into an elliptical 
galaxy. This outcome was first suggest ed by the early 
simulations of iToomre fc Toomrd ()1972l ) a nd later con- 
firmed by various other i nvestigators (e.g.. lBarnesHl988l 
iBorne fc Richstonel Il99ll ). Evidence gathered over the 
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last few decades conclusively shows that at least a frac- 
tion of elliptical galaxies were indeed for med by mergers 
(see, e.g., the review by ISchweizei1ll998t ). Such evidence 
includes the finding from ground-based imaging that on 
global scales many merger remnants are well fit by the 
i? 1 / 4 luminosity profiles that are characteristic of ellipti- 
cal galaxies (e.g.. ISchweizerl ll982|; I ffibbard fc Yur]|1999l: 
iRothberg fe Josephll2004l l 

Previous studies of spiral-spiral mergers have 
focused extensively on the tidal tails (e.g ., 
iHibbard fe van Gorkoml 119961: iHibbard et all [2001), 



globular clusters ( e -K-i IWhitmore et alj 1 1993 : 
I Whitmore fe Schweizerl Il995t ISchweizer et al.1 11996 . 



2004), young star clusters (e.g., IWhitmore et al.l 



1999, 



120051) and dwarf g alaxies that may h a ve formed during 
the m erger (e.g., iBraine et afl 120011 : IWeilbacher et all 
|2003| ). Less attention has been paid t o the nuclei, 
howe ver (although see ISchweizeil 119961 : IScoville et al.l 
l2000f ). This is in contrast with th e case for, e.g., the 
nuclear regions o f spiral galaxies (ICarollo et al.l 1997: 



Boker et al.ll2002D and ellipticals (e.g. JLauer et al.1 



1995; 



Ravindranath et all 120011; iRest et aim OOll; lLaine et all 



2003al iFerrarese et al.l l2006t l . which have been studied 



at high resolution using the Hubble Space Telescope 
(HST). In fact, one of the most important physical 
processes that occurs during the merging process is the 
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inflow of gas into th e nucle a r regio n (e.g.. [Noguchilll988t 



Mihos fc Hernquistl 11994 119961 : iBarnes fc Hernquistl 
1996|) . To study the question of whether and how merg- 



ers convert the structure of the nuclear regions from 
disk galaxy bulges or pseudo-bulges into the density 
profiles seen in the centers of elliptical galaxies, it is 
therefore important to investigate the nuclear regions at 
t he highes t achie vable angular resolution. 

iToomrd (|1977f) assembled a small sample of prospec- 
tive merger systems in various stages of merging that 
later became known as the Toomre sequence. We in- 
vestigated various properties of this optically selected 
merger sample with 77ST/WFPC2, in parti cular the re- 
sponse of the nuclei to the merging process (|Laine et al.l 
2003b, hereafter Paper I). However, in the optical sev- 
eral of the nuclear regions are hidden behind dust, which 
makes identification of the nuclei difficult, if not impos- 
sible. Near-infrared (NIR) observations are much less 
affected by dust, and thus are better suited to identify- 
ing the location of the nuclei within the often complex 
circumnuclear environments of the host galaxies. 

Previous investigations have studied a r ange of basic 
properties of various merger samples (e.g.. I Joseph et alJ 
1984 lLonsdale et all 11984 IStanford fc Bushousd 119911: 



Bushouse fc Stanford! \199% iRothberg fe Joseph! 12004 



Most of the Toomre sequence galaxies have been stud- 
ied previously from the ground at NIR wavelengths: 
NGC 4 038/39 dBushouse fc Wer^ ll990l: iHibbard et all 
[20011 : iKassin et al.l I2003D : NGC 4676 (J. Hibbard 
unpu blished); NGC 7592 dBushouse fc Werner! 

1990); N GC 6621/22 dBushouse fc Werner! I1990D; 



NGC 3509 dBushouse fc Stanford! I1992D; N GC 520 
(|Bushouse fc Stanford! Il992t iKotilainen et all 1200 If ); 
NGC 3256 (jKotilainen et al.lll996f h and more recently. 



all t he four late-stage mergers (|Rothberg fc Joseph! 
2004). However, these observations cither lacked ade- 
quate spatial resolution to reveal the detailed structure 
of the nuclei, or were aimed at investigating properties 
on larger scales. In order to explore the nuclear regions 
themselves, the high spatial resolution of the HST (5-50 
pc at the distance of the Toomre sequence galaxies) is 
imperative. 

We present here a study of the entire Toomre sequence 
with HST /NICMOS. This work is part of a high spatial 
resolution multi-instrument HST survey of the nuclear 
regions in the Toomre sequence (cf. Table [I} , including 
WFPC2, NICMOS and STIS observations. Earlier (see 
Paper I) we reported on the nuclear morphologies and 
luminosity densities based on WFPC2 bro ad- and nar- 
rowb and observations. In a future paper (jRossa et al.l 
120071 hereafter Paper III) we will present stellar popula- 
tion studies based on STIS G430L spectra. In this pa- 
per we focus on the nuclear morphology and photometric 
properties as revealed by the NICMOS observations. We 
make a detailed comparison with our previous WFPC2 
observations and other existing multi-wavelength obser- 
vations, and discuss the results in the context of the evo- 
lution of merging galaxies and the formation of ellipticals. 

Our investigation is the first comprehensive high spa- 
tial resolution study of an optically selected sample 
of mergers. Mer ging galaxies often have high in- 
frared luminosities d Sanders et al.l Il988f ) , sometimes in 
excess of Lir > 10 L @1 and samples of LIRGS (Lu- 
minous InfraRed Galaxies) and ULIRGS (Ultra Lu- 



minous I nfraRed Galaxies) are d ominated by merging 
systems (|Sanders fc Mirabel Il996f ). The present opti- 
cally selected sample thus provides a useful compari- 
son to 7JST/NICMOS studies of merger remnants se- 
lected only on the basis of their infrare d luminosity (e.g., 
iScoville et aTI[200Tll: [Borne et a.l.lf200(Th . 

This paper is structured in the following way: in §[3]we 
present the observations and the data reduction proce- 
dures. In §[3] we show results on the nuclear morphology 
and color mapping for each individual system separately 
(grouped in three merger stages), and compare them with 
other multi-wavelength studies. In § [4] we present and 
discuss photometry, colors, and radial color profiles for 
the Toomre nuclei. In § [5] we present -ftf-band surface 
brightness profiles, and discuss the results in the context 
of the evolution of merging galaxies and the formation 
of early- type galaxies. In § [5] we summarize the main 
results of the paper. 

2. SAMPLE AND OBSERVATIONS 

2.1. Sample 

Our sample is the Toomre sequence of relatively nearby 
(D < 122 Mpc) strongly interacting and merging galax- 
ies, ordered according to the es timated degre e of com- 
pleteness of the merging process (|Toomrelll977t ). The ba- 
sic properties of the eleven systems are listed in Table [TJ 
Distances in the table and throughout this paper are 
based on a Hubble constant of H — 75kms _1 Mpc -1 , 
save for the distance to the Antennae which is based 
on ob servations of red giant branch stars (jSaviane et alJ 
I2004h . The systems span a wide range of dynamical 
phases, from galaxies early in the merging process (e.g., 
NGC 4038/39) to late stage systems (merger remnants; 
e.g., NGC 7252). A detailed description of the optical 
morphology of the nuclear regions is presented in Pa- 
per I. We emphasize that the Toomre sequence was cho- 
sen originally based on the visual inspection of photo- 
graphic plates, mostly thos e presented in the Atlas of Pe- 
culiar Galaxies (|Arplll96ll ). and the rough mer ger stage 



has b een derived from the degree of coalescence ( Toomre! 
Il977f) . Although the very early and very late stages can 
be distinguished quite unambiguously, an exact place- 
ment in a sequence may be possible only after deriv- 
ing the dynamical history of the merger through detailed 
simulations. 

2.2. Observations and Data Reduction 

All NIR images discussed in this paper were obtained 
with the Near Infrared Camera and Multi-Object Spec- 
trometer (NICMOS) aboard HST, which uses a 256 x 256 
pixel HgCdTe array. We used the NIC2 camera which 
has a field of view of 19'.' 2 x 19'.' 2 and a pixel scale of 
0'.'075pix _1 . We observed seven of the eleven Toomre 
systems in three filter passbands (F110W, F160W and 
F205W), which resemble but do not quite match the 
ground-based J, H, and K bands. Additionally, we 
made use of archival NICMOS images of the remaining 
four late-stage Toomre galaxies (NGC 2623, NGC 3256, 
NGC 3921 and NGC 7252), which originate from the fol- 
lowing GO programs: NGC 2623: #7219 (P.I.: N. Scov- 
ille); NGC 3256: #7251 (P.I.: J. Black); NGC 3921 and 
NGC 7252: #7268 (P.I.: R. van der Marel). All NIC- 
MOS observations were performed in the MULTIAC- 
CUM observing mode. We note that for NGC 2623 and 
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NGC 3256 no F205W observations have been carried out, 
so we used available F222M observations instead. Fur- 
thermore, for NGC 3256 no F110W observations were 
available. A complete listing of all the observations is 
given in Table [2l 

Our observations were obtained shortly after NICMOS 
was reactivated through the installation of a cryocooler 
(installed on March 7, 2002). After this reactivation, the 
NICMOS aperture locations shifted slightly with respect 
to their prior locations in the focal plane. However, the 
aperture location table used by the telescope pointing 
system was not updated until later in the year. As a re- 
sult, some of our observations showed an offset of a few 
arcseconds between the target and the intended location 
on the detector. This generally was not a problem, since 
in almost all cases the target still fell on the detector. 
Only our observations for NGC 7592 were somewhat ad- 
versely affected, as discussed in § 13.21 below. 

The data obtained in the F110W and F160W pass- 
bands were pipeline-reduced, including the basic reduc- 
tion steps such as dark and flatfield correction using the 
STSDAS 9 tasks calnica and calnicb. We used a dither 
pattern for the observations in order to allow for good 
cosmic ray removal and correction of instrumental blem- 
ishes such as bad pixels and the "grot" (areas of reduced 
sensitivity due to paint flecks that have fallen on the 
detector). For the F205W observations a slightly differ- 
ent observing strategy was used, in order to subtract the 
thermal background emission from the telescope. This 
is only an issue for observations at or longward of the 
-fT-band. We used a small scale dither pattern in combi- 
nation with a chop between the object and background 
positions. The chop positions were separated by either 
60" or 90" (depending on the extent of the galaxy) from 
the galaxy and contain little or no emission from the 
galaxy. As this combined strategy is currently not sup- 
ported in the pipeline reduction, we had to manually 
adjust a few steps, and edit the respective table-files in 
calnica, before we could go on with processing the images 
with calnicb. For further details concerning the general 
reduct ion steps we refer to the NICMOS Data Hand- 
book ({Dickinson et all [2002V Throughout the remain- 
der of the paper, we refer for simplicity to the F110W, 
F160W and F205W/F222M filterbands as the J, H, and 
-ftT-bands. 

For the photometric comparison with our WFPC2 data 
and in order to create color maps, we had to register and 
resample the images so that the WFPC2 and NIC2 im- 
ages are on the same spatial scale. In most cases we used 
only the PC chip of the WFPC2 (field of view: 35" x 
35", as compared to 19"2 x 19"2 for NIC2) on which all 
the nuclei save for the Antennae were centered. We used 
the STSDAS task reflresamp to perform the necessary 
steps, aligning all images to the NIC2 images. Rota- 
tions, translations and scale differences were taken into 
account using a combination of header information about 
the telescope pointing and orientation, and matching ob- 
ject positions in the field. We measured the positions of 
at least two reference points (e.g., stars, compact clusters 
or H II regions) in both the WFPC2 and NIC2 images 

9 The Space Telescope Science Data Analysis System is a prod- 
uct of the Space Telescope Science Institute, which is operated by 
AURA for NASA. 



and checked the final registering by blinking the images. 
In the case of NGC 520 S and NGC4676 A we note that 
due to the lack of any reference points in the NIR that 
had visible counterparts in the optical, the final register- 
ing between the WFPC2 and NICMOS images may be 
off by 2-3 pixels. We did not correct for differences in 
the FWHM of the point-spread function (PSF) between 
the optical and infrared datasets; this may cause small 
artifacts in the central ~ 0"2 of the color maps presented 
below. However, none of the main conclusions of our pa- 
per are affected by these artifacts. 

Magnitudes were obtained from count rates (CR) using 
^Filter = -2.5 log(CR) + ZP Filtcr , where ZP Filtcr is the 
magnitude zeropoint of the individual filter. Because the 
NICMOS filters do not correspond exactly to the </, H 
and if-bands, there are color terms in the zeropoints of 
the photometric transformation equations. These color 
terms can vary throughout each galaxy, due to stellar 
population gradients and the generally very inhomoge- 
neous distributions of dust. To calculate and correct for 
the color terms on a pixel by pixel basis is not straight- 
forward. We have therefore taken a simplified approach. 
First, we adopted an initial guess zeropoint that depends 
only on the filter. These zeropoints were based on the 
color terms appropriate for a 10 9 y ear old population 
(as given bv lBruzual fc Chari ot 2003) without additional 
reddening. The color terms were calculated using the 
calcphot task in the STSDAS package synphot. We then 
refined the zeropoints on a galaxy-by-galaxy basis. For 
each galaxy we calculated (as discussed in § 0]) the count 
rates in several apertures of fixed sizes. These were com- 
bined with the initial-guess zeropoints to obtain initial 
guesses for the J, H and K magnitudes. We then used 
these magnitudes to estimate the colors in the aperture, 
and subsequently we used these colors to improve the 
zeropoint estimates. 

From a variety of tests we infer that this iterative 
scheme yields magnitudes that are accurate to ~ 0.05 
mag in each of the J, H and if -bands. This is generally 
much smaller than the variations in color either within a 
given galaxy, or between different galaxies. When quot- 
ing aperture magnitudes in § [4] we use the zeropoints 
calculated specifically for each given combination of fil- 
ter, galaxy and aperture. When displaying the images, 
constructing surface brightness profiles, and constructing 
color maps, we adopt a single zeropoint for each combi- 
nation of filter and galaxy, by averaging the zeropoints 
calculated for 100 pc and lkpc apertures. Our calibrated 
magnitudes are in a system defined by the </, H and 
if-band filter transmi ssion curves of the Bessell system 
ijBessell fc BrettJ[l988h . 

3. NUCLEAR MORPHOLOGY 

3.1. Introduction 

Figures CD (NGC 4038) through Ql (NGC 7252) show 
images of the nuclear regions of the Toomre sequence 
galaxies. Most figures show three panels arranged in the 
following order: the left panel in each figure shows the 
WFPC2 F555W (hereafter U-band) image from Paper I; 
the middle panel shows the NICMOS if-band image from 
the present study; and the right panel shows the V — K 
image. A comparison of the V and K bands provides the 
longest baseline in color space available from our HST 
data. We also created color images of J — H . J — K . 
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and H — K, composed entirely of the NICMOS data. 
However, these generally show the same structure as the 
V — K images, at lower contrast, so we do not present 
them here. Also, most galaxies reveal very similar mor- 
phology in the J, H, and K band, so we only show the 
K images. 

In this section we describe the nuclear morphologies 
on a galaxy-by-galaxy basis, focusing on the informa- 
tion revealed by the J, H, and K images. Of course, 
many of the Toomre sequence galaxies have also been 
observed at other wavelengths. For example, seven 
of the eleven systems have X-ray imaging from either 
Chandra and/or XMM-Newton in the X-ray regime, 
while NGC 6621/22 was studied with ROSAT/PSPC; 
only NGC 3509, NGC 7592 and NGC7764A have not 
yet been studied in X-rays. Many of the Toomre se- 
quence galaxies have also been studied with moder- 
ate angular resolution in H I and in the radio con- 
tinuum and in CO lines (e. g., Norris & Forbes 119951 : 
Hibbard fc van Gorkoml 119961: iGe orgakakis et all [2000; 
Hibbard et all 120011 : llono et al.ll2005[ ). All these obser- 
vations at other wavelengths have provided a wealth of 
information. However, we highlight in the discussions 
below only those aspects that have a direct impact on 
the interpretation of features observed in the NICMOS 
images. 

One advantage of our high spatial resolution NICMOS 
data is that it allows better identification of the nuclei 
in the Toomre sequence galaxies than do either ground- 
based or optical data. In this context we define "nucleus" 
loosely as a light concentration that has properties con- 
sistent with it being the center of the galaxy itself, or one 
of its merging progenitor galaxies. Table Q] lists the nu- 
clei identified in our data, including their right ascension 
and declination. These positions are measured from the 
NICM OS data, and have an absolute accuracy limited to 
~ l'/6 (jPtak et al.ll2006h . due to guide star catalog un- 
certainties and uncertainties in the relative positions of 
the HST instruments and Fine Guidance Sensors (FGS) 
in the HST focal plane. For the NICMOS imaging of 
NGC 3921, we have some evidence that the astrometry is 
in fact off by several arcseconds, and believe that the ab- 
solute positions inferred from the WFPC2 data, as listed 
in Paper I, are more accurate. While we typically have 
not marked in Figures [Tl [TBI the positions of the nuclei we 
identify, their positions can always be easily located (if 
they are not already obvious from the images) by com- 
bining the positions in the table with the coordinates 
listed along the image axes. 

In order to unambiguously identify the nuclei of the 
Toomre sequence, it is important to compare the derived 
infrared positions of the nuclei with measurements at 
other wavelengths that are relatively insensitive to dust, 
such as X-rays, radio continuum and CO emission. Ta- 
ble [3] lists the offsets between the absolute positions in- 
ferred from our NICMOS data and those inferred from 
data in these other wavelengths (where available). In 
general the positions from different wavelengths agree 
within the uncertainties; in the discussions below we ad- 
dress the nuclear positions only in those cases where the 
agreement is not good. 

Many of the Toomre galaxies have quite complex nu- 
clear morphologies in the NICMOS images. Some inter- 
esting features are obvious when blinking or displaying 



images at various contrast levels on a high resolution 
computer screen, but can be much harder to identify 
and interpret on the figures reproduced here. The discus- 
sions presented below are therefore intended to guide the 
reader through the most relevant features. Descriptions 
such as "compact" , "concentrated" , and "extended" are 
qualitative in nature and refer only to the visual impres- 
sion obtained from the images. The term "circumnuclcar 
region" is used loosely to denote the region in the im- 
mediate vicinity of the nucleus. No particular scale is 
necessarily implied, but we generally use this term to de- 
note only a small fraction (say, 10%-30%) of the size of 
the field-of-view. The NIC2 field of view corresponds to 
7.5 kpc at the median distance of the sample galaxies, so 
the circumnuclear region does not generally extend be- 
yond ~ 1 kpc from the nucleus. When referring to the 
entire NIC2 field of view (which is generally much smaller 
than the full size of the sample galaxies) we refer to the 
"central region" of the galaxy. More quantitative mea- 
sures of the nuclear structure, such as detailed aperture 
photometry, color index profiles and surface brightness 
profiles of the nuclei, are presented and discussed in §§[4] 
and[5l 

3.2. Early-stage systems: NGC4038/39 (The 
Antennae), NGC 4676 (The Mice), NGC7592, 
NGC7764A 

The Antennae 10 were imaged in two separate point- 
ings, one centered on the nucleus of NGC 4038 (Fig- 
ure[TJ) and the other centered on the nucleus of NGC 4039 
(Figure [2]). As a result of the NICMOS aperture loca- 
tion offset discussed in § 12. 21 the images were not cen- 
tered exactly on the presumed nuclear positions (deter- 
mined from previous WFPC2 imagery). The nucleus of 
NGC 4039 is an easily identified concentrated source (see 
Figure [2]). By contrast, NGC 4038 does not seem to pos- 
sess a highly concentrated nuclear region; however, a 
probable nucleus can be identified in all three NIR fil- 
terbands (marked with a box in Figure [JJ . Given the 
patchy dust obscuring the view into the nucleus in the 
optical regime, optical identifications of the nuclei can 
be misleading. In the case of NGC 4038, the nucleus 
identified in Pa per I is most likely a c ircumnuclear H II 
region (see also IWhitmore et al.1 Il999h . It is very blue 
in V — K, and is likely a foreground star- forming knot 
located in front of a region of high extinction. The nu- 
cleus identified here using the NICMOS images is barely 
visible in the optical. 

A comparison with the results obtained from data at 
other wavelengths does not provide an unambiguous an- 
swer to the question of whether the nucleus of NGC 4038 
has been correctly identified. This may not be surpris- 
ing, given the intrinsically complex nature of the nuclear 
region. We have identified the clump closest to the cen- 
ter of the overall light profile as the nucleus; this po- 

10 We use a different distance to the Antennae as compared to 
Paper I. In that paper we adopted the distance of 21.6Mpc implied 
by the observed redshift of the Antennae, neglecting any potential 
peculiar velocity of the system. Instead, we adopt here the value 
of 13.8 Mpc, which was recently d erived from measure ments of the 
tip of the red giant branch (TRGB lSaviane et a0l2004l ). While this 
distance puts the Antennae closer than has been assumed in most 
previous work, the TRGB-method has been well-demonstrated and 
the inferred distance should be more accurate than the canonical 
rcdshift-distance. 
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sition is consistent with that implied by radio contin- 
uum data, if one assumes that the nucleus coincides with 
the brightest source observed at A6cm (as assumed for 
Table [3]). However, there are several bright sources in 
the nuclear area. If the radio continuum nucleus is pre- 
sumed to reside at the centroid of these sources (as do 
iNeff fc UlvesTa d ( 2000), see their Figure 5a) then the ra- 
dio continuum and NICMOS positions are not in mutual 
agreement. A more accurate alignment of the near-IR 
and radio data might be p ossible using the technique of 
IWhitmore fc Zhangl ((2002), but we do not pursue that 
here. 

The Antennae have also been studied extensive ly in 
the X-ray regime (e.g.. iFabbiano et al.lll997l 120011 and 
follow-up papers). The NICMOS position of NGC 4038's 
nucleus is c onsiste nt with source #24 in the list of 
IZezas et"ail (|2002al fbT) . and may well be the nucleus ( as 
assumed for Table [3]). However, IZezas et"al"1 (|2002al fbT) 
actually identify their source #25 as the nucleus, offset 
in declination by about 6" from the NICMOS position. 
There is little evidence in the NICMOS images to favor 
source #25 as the nucleus, and in fact source #24 has a 
much larger X-ray flux. 

We have observed NGC 4676 in two pointings, and 
clearly detect the nuclei of both galaxies. The northern 
galaxy, NGC 4676 A (Figure [3]) , shows a highly concen- 
trated infrared source, completely obscured in the opti- 
cal, which we identify as the nucleus. The second nucleus 
(NGC 4676 B) shows a concentrated nuclear region with 
a smooth appearance (Figure |4j. 

We acquired NICMOS images of NGC 7592 in one 
pointing, as both nuclei in principle fit into the field of 
view of the NIC2 aperture (see Figure [5]). However, this 
observation was affected by the NICMOS aperture loca- 
tion offset, discussed in § 12.21 which placed the second 
circumnuclear region right on the edge of the NIC2 aper- 
ture. Although the image covers most of the circumnu- 
clear region, this region shows a very complex morphol- 
ogy and there is a chance that the true eastern nucleus 
(NGC 7592 E) is located outside the field of view. In 
contrast, the western nucleus (NGC 7592 W) is clearly 
defined and compact , and was previously identified by 
iRafanelli fc Marzianil (|1992ft to be of Seyfert 2 type. 

The nuclear region of NGC 7764A is shown in Figured 
When enlarged, all three NICMOS images (J, H and K) 
reveal what appears to be a double nucleus in this system 
(see Figure [7]). If these truly are the two merging galaxy 
nuclei, it is consistent with the placement of NGC 7764A 
in the middle of the Toomre sequence. The two puta- 
tive nuclei are separated in a nearly east-west direction 
by 0'.'43, which corresponds to 260 pc at the distance 
to NGC 7764A. This is similar to the sep aration of the 
two n uclei in the merging galaxy Arp 220 (|Scoville et al.1 
[19911 . 

3.3. Mid-stage systems: NGC6621/22, NGC3509, 
NGC 520 

The nucleus of NGC 6621 is detected in the infrared, 
and has a fainter counterpart in the optical images (see 
Figure [H]). The circumnuclear region has a complex mor- 
phology, even in the if -band. In contrast, the nuclear 
region of NGC 6622 appears much less complex, showing 
only a single concentrated source which we identify as 
the nucleus (Figure [9]). 



The nuclear region of NGC 3509 looks rather inconspic- 
uous, with only a single, well-defined nucleus (FigurefTO]). 
appearing as a concentrated but clearly resolved source. 
As noted already in Paper I, NGC 3509 has possibly ex- 
perienced a minor merger in the past, as onl y one clea r 
tidal tail is visible in deep optical images (cf. lArplll966h . 
The relatively low value of the IRAS flux density ratio at 
60/im over lOO^rn (/6o / /ioo ~ 0.33) also favors a scenario 
in which star formation has been suppressed compared 
to the other mid-stage mergers in the Toomre sequence, 
which show more vigorous signs of recent or ongoing star 
formation (cf. Figures 7 and 8 of Paper I) . 

NGC 520 was also observed in two pointings. The NIR 
observations of the southern nucleus (NGC 520 S) are 
quite striking (Figure [TT|). revealing a nucleus completely 
hidden in the optical behind an enormous dust complex. 
The J-band observations show a considerable amount of 
dust present around the nuclear region. At longer wave- 
lengths (H and K) we peer deeper through the dust to- 
ward the nucleus, but even the K image is affected by 
patchy dust structure around the circumnuclear region. 
However, the nucleus is clearly visible, appearing as an 
extended disk composed of at least seven major knots. 
The brightest knot, which we identify as the nucleus, is 
slightly extended toward the north. Interestingly, the 
nuclear disk, which appears almost edge-on in the E-W 
direction, is slightly skewed with respect to the large scale 
dust l ane. Such a mi salignment is also seen in Centau- 
rus A (jSchreier et al.lll998l ). 

Previous sub-arcsecond resolution ground-based NIR 
im aging also revealed an e xtended disk around NGC 520 
S (|Kotilainen et al.ll2001f h However, considerable sub- 
structure was only seen in Br7 observations, not in 
their if -band ima ge. Radio observations at A 2 cm 
(|Carral et al.lll99"ot) also showed an extended disk with 
five prominent knots, very similar in morphology to our 
five brightest components. More recently, high spatial 
resolution radio continuum observations at A 21.4 cm ob- 
tained with MERLIN revealed an extended disk w ith at 
least nine bright components (jBeswick et al.l [2003). also 
consistent with the NIR morphology. The much less ob- 
scured northern nucleus (NGC 520 N) reveals a bright 
centrally confined, compact morphology, embedded in an 
amorphous disk that is highly inclined to the line-of-sight 
(see Figure [T2")) . 

3.4. Late-stage systems: NGC 2623, NGC 3256, 
NGC 3921, NGC 7252 

The nuclear region of NGC 2623 shows a highly con- 
centrated light profile (Figure fT3|) . Consistent with our 
previous findings (see Paper I), only one nucleus is iden- 
tified, appearing slightly elongated in the E-W direction. 
There is no clear nucleus discernible in the optical image. 

Our analysis of the late-stage merger galaxy NGC 3256 
is based on archival NICMOS images in H and F222M 
(hereafter referred to as K; see Figure [T4"|) ; no J band 
data exists. The nuclear region (NGC 3256 N) appears 
as a concentrated and smooth s ource. Radio imaging a t 
A 3 cm and A 6 cm, performed bv lNorris fc Forbes! (pm 
identified a second nucleus 5" south of the primary 
(NGC 3256 S). Studies of molecular gas e mission also hint 
at a s eco nd nucleus a t the sa me location (jSakamoto et al.l 
[200l . iBoker et all (fl997h also imaged NGC 3256 in 
the iV-band and detected a source at this same posi- 
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tion. Based on r adio and X-ray power and their ratios 
iNeff et al.l ((2003) argue that both nuclei may be low- 
luminosity active galactic nuclei (LLAGNs). In our NIC- 
MOS data there is a compact source at the position of 
the secondary nucleus, but with fainter intensity than 
the primary nucleus. The source appears projected onto 
the middle of a spiral arm, lying in a highly reddened 
region corresponding to a local maximum in the V — K 
image. It has an extinction-corrected X-band magnitude 
of —18.2, similar to the nuclei in the other Toomre se- 
quence galaxies (see £|4.3[) . All this supports the view 
that there is in fact a highly extincted infrared counter- 
part to the secondary nucleus that was detected in radio, 
MIR, and X-ray data. 

We also analyzed archival images of NGC 3921. The 
nuclear region shows a slightly elongated core in the NNE 
to SSW direction (see Figure [T5]) . with only one nucleus 
det ected. Our findi ngs are consistent with earlier stud- 
ies (|Schweizerlll996D and with the fact that NGC 3921 is 
generally regarded as a late-stage merger remnant with 
properties re sembling elliptical galaxies. In a recent X- 
ray study bv lNolan et a l. (2004), however, the X-ray lu- 
minosity of the hot diffuse gas was found to be two orders 
of magnitude less than that typical of elliptical galaxies. 

Our images of the latest-stage merger remnant of the 
Toomre sequence, NGC 7252, show a nuclear spiral (see 
Figure tTHJ) , as already detected in the optical and de- 
scribed in our previous study (Paper I). Only one nu- 
cleus is identified in all three filter passban ds, consis- 
tent with previous ground -based studies (e.g.. ISchweizerl 
Il982t iHibbard et alj I1994D and with the general notion 
of NGC 7252 being a prototype mer ger remnant, wi th an 
estimated age of about 1 Gyr (e.g.. ISchweizerl fl998h . 

4. PHOTOMETRY AND COLORS 

4.1. Photometry 

We measured the J, H, and K magnitudes for all iden- 
tified Toomre sequence nuclei within apertures 100 pc 
and 1 kpc in size. Depending on the distance to each 
galaxy, the 100 pc aperture corresponds to 2-20 pixels 
and the lkpc aperture corresponds to 23-200 pixels. 
Fluxes were measured with the IRAF 11 task phot, cen- 
tered on the identified nuclei; the resulting magnitudes 
are listed in Table [4] and Table [5] The photometric ran- 
dom errors on the magnitudes and colors are very small 
(with a median of 0.01 mag) and are not listed in the 
tables. The actual errors are dominated by systematic 
calibration uncertainties (< 0.05 mag in each band, see 

4.2. Near-Infrared Colors 

Figure [17] shows the J — H vs. H — K color-color 
diagram for the 100 pc aperture measurements of the 
Toomre sequence nuclei. Separate symbols are used for 
the early-stage, mid-stage, and late-stage mergers, as 
indicated in Table [TJ For a direct comparison of the 
late-stage mergers with elliptical galaxies, we also show 
in the figure the sample of elliptical galaxies studied 
by ISilva fe Bothunl (|1998l ). We transformed their data, 

11 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agreement 
with the National Science Foundation. 



which were calibrated to the CIT photometric system, 
to the Bes sell system in which o ur NICMOS data were 
calibrated (|Bessell fc Brettlll988t ). We also used the re- 
spective J — H and H — K values from the GALAXEV 
database to generate and over plot the evolutionary trac k 
of an aging simple population (jBruzual fc Cha riot 2003) . 
We used the Bruzual-Charlot models with the Padua 
1994 stellar evolution library for a standard Chabrier 
initial mass fu nction (IMF) from 0.1 M to 100 M Q 
(|Chabrierll2003l ) and solar metallicity. Two of the three 
late-stage merger remnants (those which had J, H, and 
K measurements) fall clos e to the positions of t he el- 
liptical galaxies studied by [SiIiOEBothu3 ((1991) . All 
Toomre sequence nuclei are closely aligned along the red- 
dening vector direction, as expected for dusty stellar pop- 
ulations. 

Figure [TH] shows the same color-color data for the 
Toomre nuclei as Figure [T7l but focuses on a smaller re- 
gion of the color-color diagram. W e compare our data to 
groun d-based measurements from iBushouse fc Stanford! 
(1992) (also transformed to the Bessell system) for a 
larger sam ple of inter acting galaxies selected from the 
Arp atlas (|Arpl 1 1966( 1. The average and median col- 
ors are slightly redder for the Toom r e seq uence nu- 
clei than for the IBushouse fc Stanford! (|1992f ) galaxies, 
but this result is only marginally significant given the 
large scatter in both distributions. Moreover, a di- 
rect co mparison between the samp les is not straightfor- 
ward. IBushouse fc Stanford! (|1992l ) performed photom- 
etry in a fixed angular aperture of 5'.' 4, which is gen- 
erally much larger than the apertures of a fixed physi- 
cal size that we have considered for the NICMOS data. 
Also, Toomre's criteria for characterizing a galaxy as 
in teracting or merging were m ore stringent than those 
of IBushouse fc Stanfordl {1992). In any case, Figure [T5] 
shows that the majority of the Toomre sequence nuclei 
fall in the region of the J — H v ersus H — K diagram 
that is populated by the sample of IBushouse fc Stanfordl 
(1992). The Toomre sequence nuclei therefore have col- 
ors that are comparable to those seen in other interacting 
and merging galaxies. 

In principle it should be possible to constrain the stellar 
populations in the nuclei of the Toomre sequence galax- 
ies using the multi-band V, I, J, H and K photometry 
presented here and in Paper I. Such analyses have al- 
ready been perform ed for other merging galaxies (e.g., 
iPasquali et al.1 120041 ). However, in the presence of sig- 
nificant dust, as is the case here, the results are quite 
uncertain. We therefore postpone a discussion of the 
stellar populations to a subsequent paper (Paper III), in 
which we will present HST/STIS spectra obtained in the 
context of our investigation. These spectra can be used 
to investigate the populations of the nuclei with moder- 
ate obscuration. By contrast, near-infrared spectroscopy 
would be required to study the populations of the most 
heavily obscured nuclei. 

4.3. Dependence on Merger Stage 

In Figures [19] and [20] we examine color trends along 
the proposed merger sequence. The NIR J — K colors 
in apertures of 100 pc and 1 kpc are shown as a function 
of merger stage in Figure [TH The J — K colors vary 
between 1.0 and 2.1 for all but the two most heavily 
obscured nuclei (NGC 2623 and NGC 520 S), which have 
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much redder values within the 100 pc and 1 kpc apertures 
(see Tables E] and respectively). We have also plotted 
V — K colors (listed in Tables 2] and O as a function 
of merger stage (see Figure |20|) . therefore extending the 
spectral baseline to optical wavelengths. By and large, 
these yield results consistent with those seen for J — K. 
At the latest stages, the NIR colors tend to be less red 
and closer to those of normal E/SO's. Out of 15 nuclei, 
the three bluest nuclei in a 100 pc aperture are all found 
in the five latest galaxies of the sequence. While this is 
suggestive of a trend for the nuclei to become bluer with 
increasing merger stage, it is clear that there is large 
scatter in color along the sequence, and that any trend 
is weak at best and not monotonic. 

Unfortunately, the near-IR colors measured here do 
not provide much insight into variations in age along 
the Toomre sequence, or variations in age between the 
Toomre sequence galaxies and elliptical galaxies. The 
J — K color of a stellar population does not vary much 
as a function of a ge. For exampl e, a solar-metallicity 
population with a IChabrierl (|2003f) IMF has J - K = 
0.80, 0.57, 0.82, and 0.89, at ages of 10 7 , 10 s , 10 9 , 
and 10 10 years, respectively (using the Padua 1994 stel- 
lar libra ry and an IMF that ext ends from 0.1 M© to 
100 M^; iBruzual fc Chariot! I2003D . Moreover, realistic 
stellar populations for the galaxies in the present sample 
cannot produce J — K colors that are much redder than 
~ 1.0. Therefore, the redder J — K colors seen in Fig- 
ure [19] must be due to dust extinction. As a corollary, 
the trend that the early-stage mergers are redder than 
the late-stage mergers merely indicates that the former 
are dustier than the latter. 

Models of merging galaxies with black holes often 
invoke feedback from an AGN to disp erse obscuring 
dust and shut off star formation (e.g., Hopkins et al. 
2005). Combined with the hierarchical formation of 
structure over a Hubble time, such models can suc- 
cessfully explain a variety of observed properties re- 
lated to galaxy for mation and quasar evolution (e.g., 
iHopkins et"alll2006f ). The existence of AGN-driven feed- 
back is supported by observations of individual galaxies 
such as NGC 6240, which has two AGN and a stro ng su- 
perwind (jKomossa et al.ll2003l : lGerssen et al.ll2004f ). Our 
finding that the dust content of the Toomre sequence nu- 
clei has a decreasing trend with merger stage is qualita- 
tively consistent with the predictions of feedback models. 
While the Toomre sequence does not currently host many 
strong AGN (see Table [lj , this is not difficult to recon- 
cile with the feedback models if bright A GN stages are 
eithe r very short-lived or highly obscured (jHopkins et al.l 
2005). It will be valuable to use numerical models to 
make quantitative predictions for the dust and gas con- 
tent of interacting and merging galaxies as a function 
of merger (st)age, and compare these to results such as 
those obtained here. 

As in our previous WFPC2 study, we also calculate the 
total luminosity within the 100 pc and 1 kpc apertures for 
the sample galaxies, in this case for the if -band. The K- 
band luminosity {Lk) is shown in Figure [2"T1 as a function 
of merger stage, where the luminosities have been cor- 
rected for dust extinction. To perform this correction, 
we first assume that the stellar population has an intrin- 
sic J — K m 0.85. The Jf-band extinction for each nu- 
cleus can then be estimated from Tables|4]and[5]using the 



equation A k = 0.659(J — K — 85), based on the extinc- 
tion law of iRieke &; Lebofskvl ([l985). 12 Our results are 
not significantly affected for other realistic choices for the 
intrinsic J — K color. Figure [2"T1 shows that Lk (within 
an aperture of fixed physical size) tends to increase with 
advancing merger stage. To statistically test this trend, 
we calculate the Spearman rank-order correlation coeffi- 
cient (r s ) for each of the lOOpc and lkpc datasets. In 
both cases we find that r s = 0.61, indicating a 2a corre- 
lation (99% CL). We note that the Antennae drive much 
of this correlation; removing the Antennae nuclei from 
the analysis lowers the correlation coefficient to r s = 0.5, 
indicating a 1.5a correlation (93% CL). We caution, how- 
ever, that this analysis assumes that there are no sys- 
tematic errors which depend on the merger stage. While 
not shown here, we have done the same analysis for the 
fractional nuclear luminosity (i.e., Lk,uucI 'Lx,tot) as a 
function of merger stage and find a similar trend. Be- 
cause we have corrected for extinction, the underlying 
cause of this trend must be either an increasing stellar 
density with advancing merger stage, or a younger av- 
erage stellar age with advancing merger stage (or both). 
We cannot discriminate between these two possibilities 
on the basis of the near-IR imaging data alone. Avail- 
able far-infrared data also cannot address this issue, since 
the global star formation rates measured b y IRAS (e.g., 
I.Tosenh fc Wright! fl985l iCasoli et al.lfl99lh ' are not nec- 
essarily a good indicator of the star formation rates in 
the nuclei (Kcnnicu tt et al.lll987f K As a caveat we would 
like to note that the observed trend of increasing Lk as 
a function of advancing merg er stage is base d on the ini- 
tially assumed merger stage (|Toomrelll977|) . which may 
not necessarily reflect the actual merger (st)age for all 
of the eleven merger systems. Nonetheless, we use the 
pre-defined merger stage (which is a fixed value) rather 
than the merger age. While the very early- and very 
late-stages are well defined, an actual ranking (i.e. an 
exact placement within this sequence) may be possible 
only after deriving the dynamical history of the merger 
through detailed simulations. 

4.4. Radial Color Gradients 

The two-dimensional spatial variations of the V — K 
color in the central regions of Toomre Sequence galaxies 
are quite complex, as evident from the gray-scale rep- 
resentations in Figures fTT- HQt To obtain more easily in- 
terpretable quantitative insight we have therefore deter- 
mined the radial V — K color gradients for all of the 
nuclei, measured in concentric circular annuli centered 
around each nucleus and shown in Figure 1221 The an- 
nuli are contiguous, and grow in radius by a factor of 
1.1 at each step. The color profiles are followed out to 
a distance of 5" from each nucleus (which corresponds 
to between 300 pc and 3.0 kpc for the distances of the 
sample galaxies). No correction was made for the PSF 
differences between the V and if-bands, so V — K fea- 
tures in the central ~ 0'.'2 should not be trusted. 

The last two merger remnants in the sequence 
(NGC 3921 and NGC 7252) and NGC 3509 all have rather 
shallow color gradients, with a color change of less than 

12 For NGC 3256 we used a similar formalism for calculating Ax 
based on the H — K color, because no J-band data are available 
for this galaxy. 
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one magnitude within 5". Most nuclei have steeper gra- 
dients, with a change in color of up to six magnitudes 
within 5". The nuclei with the strongest color gradi- 
ent are mostly the mid-stage mergers. The majority of 
the nuclei have nuclear colors of — K ^ A — 5, except 
for the very dusty systems (NGC520S, NGC2623 and 
NGC 3256 S) which have much redder nuclear colors of 
V-K ~7-8 (cf. Fig. HI|). 

Several of the nuclei become considerably redder in- 
ward of about 5" from the nucleus (e.g., NGC 4676 A, 
NGC7764A, NGC 520 S, NGC 2623, NGC 3256 S), prob- 
ably indicative of a strong concentration of dust in the 
central regions. This dust complicates the identifica- 
tion of potential population gradients from the observed 
color gradients. NGC 6621 becomes bluer in its central 
~ fX'6, but the morphology of the V — K image (see 
Figure [SJ) suggests that this may be due to an absence 
of dust rather than the presence of young stars. A few 
nuclei are bluer in their central ~ 0'.'2 (e.g., NGC 520N, 
NGC 7252), but this may be due to PSF differences be- 
tween V and K rather than differences in stellar popu- 
lations. On the whole, Figure [22l certainly shows no evi- 
dence for the bluer nuclear colors that might be expected 
if the galaxies all had recent nuclear star formation. On 
the other hand, this scenario is certainly not ruled out by 
the data, given that large quantities of dust affect most 
of the color gradients. 

5. SURFACE BRIGHTNESS PROFILES 

5.1. Analysis 

Surface brightness profiles for the Toomre nuclei were 
derived primarily from the K band imaging, as they 
are least affected by dust and star formation. The im- 
ages were run t hrough 20 iterations of Lucy- Richardson 
deconvolution (|Richardsonl Il972t iLucvl Il974f) using the 
VISTA software package (|Lauerlll986l ). where the use of 
a higher number of iterations did not result in a signif- 
icantly different dec onvolved image. We used TinyTim 
(|Krist fc Hookll200ll ) to generate a PSF for each nucleus, 
modeled at the specific position of each nucleus on the 
detector array and using a K-type star for the PSF spec- 
tral d istribution. For NG C 2623 we used the drizzled im- 
ages ([Scoville et al.ll2000l ) with half the pixel size of the 
original data. In each case, the PSF was generated out 
to a radius of 3 arcseconds, which incorporates the low- 
level emission in the extended wings, and was tapered at 
the edges with a Gaussian taper. A few of the galaxies 
have nuclei that are compact enough to produce a strong 
Airy ring and diffraction spike pattern. We found that 
these artifacts were in general successfully mitigated by 
our deconvolution procedure, although some residuals do 
remain (especially for the most compact sources, such as 
NGC 2623). 

The nuclear positions were determined by finding the 
centroid in a small box, typically with a size of about 
7x7 pixels. We then fitted ellipses to the isophotes of 
the two-dimensional convolved and deconvolved images. 
The surface brightness profile was derived only for the 
nuclear region (usually out to a radius of a few arcsec- 
onds). We did not attempt to mask out dust or stel- 
lar clusters within the area where the surface brightness 
profile was derived, since the effects of dust are subtle in 
most cases and there were usually no strong point sources 
within the few arcsecond radius that could have inter- 



fered with the profile. In cases where such sources exist 
(e.g., NGC 6621), the profile was only extracted within 
the innermost nuclear area. The effects of subtle dust 
features and circumnuclear spiral arms are still visible 
as small-scale wiggles in some of the luminosity profiles 
(e.g. , NGC 3256 and NGC 7252) . 

Follo wing closely the procedure used by lLaine et"aTI 
(|2003al) . for ra dii smaller th an 0'/5 we used the VISTA 
profile task (jLauerl IT985) on the deconvolved image 
to find the surface brightness profile. This task keeps 
the center fixed and fits ellipses by sampling the light 
profile in a circle wit h a radius of 1, 2, 3, etc. pix- 
els. The snuc task (jLauerl [198(1 in VISTA was used 
on the deconvolved image between the radii of 0'.'5 and 
1" to obtain the surface brightness profile. Outside of 1" 
the effects of the PSF are negligible (although diffraction 
spikes from a very compact nucleus may still be visible 
there). At these radii we therefore used the snuc task 
on the surface brightness profile inferred from the orig- 
inal (not PSF-deconvolved) images, motivated by S/N 
considerations. The transition at 0'/5 was based on ex- 
perience which indicated that the profile task is more 
robust in the interpolation of pixels at very small radii 
(Tod Lauer, priv. comm.). The results of the analysis 
are shown in Figure [23] for twelve Toomre sequence nu- 
clei. The remaining four nuclei (NGC 4038, NGC 7592 
E, NGC7764A, NGC 520 S) were too complex or other- 
wise unsuitable for ellipse fitting and profile extraction. 
The first panel shows for comparison the profile of the 
NICMOS PSF itself, indicated by star symbols. Fig- 
ure mi also shows for each galaxy the ellipticity and po- 
sition angle of the fitted ellipses, which can potentially 
be used to trace structures such as nuclear bars and spi- 
ral arms. The random uncertainties for all profiles are 
usually smaller than the size of the plotted symbols and 
hence are not shown. 

A common analytic function to fit to the luminosity 
surface density S of elliptical galaxies is the so-calle d 
"Nukerlaw" (e.g.. lLauer eFal1fl995l: iBvun et al.lll99l : 

E(r) = S 6 2^ {r/r h y~< (1 + [r/r b ] Q )^. (1) 
The surface brightness fi follows from 

/i[mag/n"] = -2.5 log E(L Q /pc 2 ) + 21.572 + M Q , (2) 

where the sola r absolute magnitude in t he -ftf-band is 
M k ,q = 3.28 (jBinnev k. Merrifieldl fl998l ) . The Nuker 
law represents a power-law with a break at the radius 
rb, with the parameter a measuring the sharpness of the 
break. The asymptotic power-law slope is 7 at small radii 
and (3 at large radii, while is the luminosity surface 
density at . The best-fitting Nuker laws are plotted in 
Figurel2"3"las solid curves, and the fit parameters are listed 
in Table O where is derived from the fitted parameter 
S 6 . 

As a diagnostic in the discussion below we generally 
use the quantity To. 15, defined as the power-law slope 
at a fixed, resolved radius of 0'.'15. This is a slightly 
more robust measure of the central brightness profile 
slope than the fit parameter 7, although the t wo are 
generally in good agreement (|Laine et al.ll2003al ). Both 
quantities are listed in Table [51 Galaxies with slope 
> 0.5 are traditionally called power-law galaxies and 
thos e with slope 5- 0.3 are traditionally called core galax- 
ies (|Faber et alj 119971) . Whether these values bracket 
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truly distinct classes of galaxies or whether there is sim- 
ply a continuum i n cusp slopes among galaxies contin- 
ues to be debated (iRavindranath et al JI200 It iRest et al.l 
l2001tlLauer et al.ll2005HFerrarese et al.ll2006D . Here, we 
merely adopt the accepted terminology, without neces- 
sarily implying that galaxies come in two physically dis- 
tinct classes. 

As a consistency check we also derive surface brightness 
profiles for a few galaxies from the J- and ii-band im- 
ages, yielding very similar results to the iT-band (apart 
from an overall color difference). As the J- and ii-band 
data have a somewhat narrower PSF than the iT-band 
data, the good agreement therefore indicates that un- 
certainties in the iT-band surface brightness profiles due 
to residual problems associated with PSF deconvolution 
are small. In the remainder of the paper we discuss only 
the surface brightness profiles derived from the iT-band 
images. 

5.2. Individual Galaxy Profiles 

The profiles for the nuclei NGC3256N, NGC3509, 
NGC 4039, NGC 4676 B, NGC 7252 and NGC 7592 W are 
essentially pure power-laws with an almost constant 
slope through the fitted region. The r .i5 values vary 
between 0.94 and 1.58, therefore belonging comfortably 
to the regime of power-law galaxies. Some deviations 
from the best fit are seen in the profiles of NGC 3256, 
NGC 7252, and NGC 7592, most likely due to intrinsic 
structure in these nuclei. The first two of these galaxies 
have circumnuclear spiral structure, whereas NGC 7592 
has an end-brightened bar-like structure. In all cases, 
residual structure from the diffraction pattern may play 
some role as well. Either way, in all three cases the over- 
all slope remains relatively unchanged from large to small 
radii within the fitting area. 

The r .i 5 values of NGC 2623, NGC 3921, NGC 4676 
A, and NGC 6621 are large enough to place them in the 
realm of power-law galaxies. However, the surface bright- 
ness profiles of these galaxies do show more structure 
than a pure power law. NGC 2623, for example, shows 
several wiggles in its profile. Nonetheless, the bright- 
ness profile is steep throughout the entire nuclear region, 
with To. 15 = 1.42 and no break toward a shallow core. 
The profile of NGC 3921 appears to have a shallow break 
around 0'.'25, but with an inner slope (ro.15 = 1.13) that 
is still in the power-law regime. NGC 4676 A has a clearer 
break in the power-law slope, around 0'.'43, but its To. 15 
value of 0.65 also still places it comfortably among power- 
law galaxies. Some curvature around 0'.'4 is seen in the 
profile of NGC 6621, but the To. 15 value of 0.68 indicates 
that this nucleus too belongs to the power-law group. 

NGC 6622 is the only galaxy for which the profile turns 
over to a power-law slope 7 = 0.24 that is characteris- 
tic of core galaxies. However, the turn-over happens in- 
side of ff.'l, the characteristic resolution of our data. So 
while this galaxy may have a core-type profile, it is cer- 
tainly not an unambiguous case. In the following we do 
not discuss the parameter of the NUKER profile fits. 
This parameter generally has physical meaning only for 
galaxies with well resolved cores, and none of the Toomrc 
sequence nuclei fall in this category. 

The only nucleus that does not fit within the class 
of either the core or the power-law galaxies is one nu- 
cleus in the mid-stage merger NGC 520 N. This nucleus 



has a very shallow profile out to several arcsec from the 
center, but with a strong upturn in the central ~ 0'.'15. 
The bright central source responsible for this upturn is 
clearly visible in the if -band image (Figure [P2f . This 
profile and nuclear morphology are very similar to what 
is see n in the centers o f many late-type spiral galaxies 
(e.g.. IBoker et al.l l2002). where it has been established 
that the upt urn is due to the presence of a nucle ar star 
cluster (e.g.. IBoker et alJl2004HRossa et al.ll2006l ). Since 
more than half of the spiral galaxies have such nuclear 
star clusters, it would not be surprising to find one among 
the Toomre sequence galaxies as well. 

Because the galaxies at the earliest stage in the Toomre 
sequence have not yet experienced the full merger evolu- 
tion, their nuclei would be expected to still have bright- 
ness profiles similar to those of spiral galaxies. At 
HST resolution, the "classical" i? 1 / 4 bulges in early-type 
spirals tend to have power -law brightness pr ofiles with 
slopes 0.5 < 7 < 1 (e.g., ISeigar et al.l I2002D . This is 
indeed similar to what we find in the earliest stage merg- 
ers of the Toomre sequence (NGC 4039, NGC 4676 A, 
NGC 4676 B). 

5.3. Comparison to Elliptical Galaxies 

As merger remnants may evolve into elliptical galax- 
ies, it is of interest to compare the nuclear properties 
of these two classes o f galax ies. We used the results 
of iRavindranath et alj (|2001l hereafter R01) as a suit- 
able comparison sample of elliptical and SO galaxies. 
While many other HST studies have addressed the sur- 
face brightness profiles of early-type galaxies in the opti- 
cal regime, the R01 study is one of the few that use ob- 
servations with NICMOS in the near-infrared. We trans- 
formed their if -band results to the iT-band assuming a 
fixed color H — K = 0.215, as appropriate for normal 
early- type galaxies (see Figure fT7|) . 

Previous studies have found that the properties of 
the nuclear brightness profiles of elliptical galaxies cor- 
relate strongly w ith the total galaxy luminosity (e.g., 
iFaber et al. 1997). To compare the properties of merger 
remnants and elliptical galaxies we therefore calculated 
the total absolute iT-band magnitudes (Mk) for both 
our own sample and that of R01, using the apparent -fT- 
band magnitudes from the Two Micron All Sky Survey 
(2MASS) and distances based on the observed systemic 
velocity. The results for our own sample are listed in 
Table [I 

Figurel2~4"lshows To. 15, the power-law slope at 0'.'15, ver- 
sus Mk- Optical studies have found that galaxies with 
M v > -20.5 (M K > -23. 6) 13 are generally power-law 
galaxies, that galaxies with M v < -22.0 (M K < -25.1) 
are generally core galaxies, and that galaxie s of interme- 
diate luminosities can be of either type (e.g.. lFaber et al.l 
1997). The data from R01, shown in Figure l24l support 
these trends (e.g., beyond Mk ~ —24.7 there are only 
core galaxies in their sample, but no power-laws). While 
these trends appea r more pronounced in some studies of 
other samples (e.g.. lFaber et al.l fl997) . it is clear that the 
correlation has significant scatter. 

Figure [2H also shows the Toomre sequence nuclei in 

13 The absolute A"-band magnitudes listed in parentheses as- 
sume V — K = 3.1, a s typical for an old stellar population 
HBruzual k. CharloUl2003l) . 
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comparison to the E/SO sample. As discussed already in 
§ 15. 2\ the Toomre nuclei correspond almost exclusively 
to power-law profiles; only NGC6622 may possibly be 
a core-type galaxy. The Toomre sequence galaxies oc- 
cupy the range of absolute magnitudes Mk in which 
both power-law E/SO and core-type E/SO galaxies can 
be found; it is therefore no surprise to find power-law 
profiles among the nuclei of the Toomre sequence galax- 
ies. However, it is of interest that so few of the nuclei 
have a core- type profile and that many have To. 15 values 
that are high relative to E/SOs. We comment on this be- 
low in the context of popular theories for the formation 
of cores in E/SO galaxies. 

Additionally, we have studied the luminosity surface 
density £ at 100 pc as a function of the total absolute 
.fT-band magnitude, shown in Figure 1251 The majority 
of the Toomre nuclei have a higher luminosity surface 
density than the E/SO's, although there are a few in the 
overlapping region. The most luminous Toomre galaxies 
also tend to have the highest surface luminosity density. 
Note that neither sample was corrected for dust extinc- 
tion in the if-band (see § 14.31 for a discussion of the typ- 
ical amount of extinction for the Toomre nuclei). How- 
ever, Figure [T71 shows that the Toomre nuclei have more 
dust extinction than typical E/SO galaxies. Therefore, 
any correction for dust extinction would in fact increase 
the systematic difference in luminosity surface density 
between the two samples. 

Figure [24] shows that, on the whole, the Toomre Se- 
quence nuclei tend to have slopes that are significantly 
steeper than what is seen in E/SO galaxies. The four 
latest-stage merger remnants all have 1.1 < To. 15 < 
1.4. In principle this could be a resolution effect, as 
the median distance of the Toomre galaxies is 77.8 Mpc, 
versus 17.0 Mpc for the galaxies in the sample of R01. 
Early-type galaxies in general tend to have steeper log- 
arithmic surface brightness profiles at larger radii, so 
that at a fixed angular scale, one would naturally ex- 
pect to find somewhat steeper slopes in a more distant 
sample. However, Figure [25] shows that this is unlikely 
to be the whole explanation in the present case. The 
Toomre nuclei tend to have a higher luminosity surface 
density than the E/SO's at fixed physical radius (i.e. at 
100 pc), and thus must be intrinsically brighter at small 
radii. It then seems plausible that they have intrinsi- 
cally steeper surface brightness slopes as well. A similar 
result was re cently repor t ed in an Z/5T/WFPC2 opti- 
cal study bv lYang etaTI (|2004h . who found that E+A 
galaxies (shown to be evolved merger remnants) had 
power-l aw profiles a nd hi gher surface brightnesses than 
E/SOs. lYang et all (|2004l ) attributed the higher surface 
brightnesses to the recent star formation in E+As and ar- 
gued that, once the A stellar population faded, the E+A 
galaxies would have surface brightness profiles typical of 
E/SO s. Our results also connect naturally with the find- 
ing of lRothberg fc Josep h (2004) that many merger rem- 
nants have higher luminosity densities near their centers 
than would be expected on the basis of inward extrapo- 
lation of the surface brightness profile at large radii. 

The steep if-band profiles and high luminosity sur- 
face densities seen in the Toomre galaxies cannot be 
attributed to the presence of bright non-thermal point 
sources. We measure the profile slope at 0'.'15 and the 
luminosity surface density at 100 pc; the latter corre- 



sponds to 0'.'26 at the median distance of the sample 
galaxies. Therefore, both of these quantities are de- 
fined outside of the direct influence of a central non- 
thermal point source, if one existed. This is true in 
particular because we measure these quantities on the 
deconvolved profiles, not the observed profiles. More- 
over, there is little independent reason to suspect the 
presence of non-thermal point sources. Most of the sur- 
face brightness profiles do not show an upturn near O'.'l. 
The only exception is NGC520N, but for this galaxy a 
nuclear star cluster is the likely cause. While some of 
the galaxy images (e.g., NGC2623 and NGC3921) do 
show the NICMOS PSF diffraction pattern around the 
galaxy center, this merely implies that the galaxy light 
must be very centrally concentrated (as indeed our anal- 
ysis confirms), not that a central point source need be 
present. Finally, non-thermal point sources at optical 
wavelengths are generally associated with powerful Ac- 
tive Ga lactic Nuclei (AGN) and r adio galaxies in par- 
ticular (jVerdoes Kleijn et al.l l200"2] ) . By contrast, most 
of the Toomre galaxies do not host strong AGN (see Ta- 
ble [J). 

The results of our surface brightness profile analysis 
have a natural explanation in the context of our un- 
derstanding of the merger process. As two galaxies in- 
teract, gravitational torques remove ang ular momentum 
from the gas and drive n u clear inflows (iNoguch: 119881 : 
iMihos fe Hernquistl Il996t iBarnes fe Hernquistl |19"96¥ . 
This transports gas to the centers of the interact- 
ing galaxies, and subsequently to the center of the 
merger remnant, triggering nuclear starbursts. The 
newly formed stars increase the stellar density in the 
nucleus and steepen the surface mass density profile 
(|Mihos fe Hernquistlll994l ). As the young stellar popula- 
tion will have a lower mass-to-light ratio than the older 
stars in the galaxy, the steepening of the projected lu- 
minosity profile will be even stronger. All this is quali- 
tatively consistent with the results shown in Figures |2~41 
and [25] The formation of new stars in a progressing 
merger may also be responsible for the observed trend 
of increasing nuclear luminosity with advancing merger 
stage shown in Figure [21] NGC 3921 and NGC 7252 are 
known to have "post-starburst" A-star features in spec- 
tra of their nuclei, which provides dire ct evidence that 
stars recently formed in these mergers (|Schweizerl 119821: 
iLiu fe Kennicuttlll995a llb1). 

To quantitatively compare the mass densities of ongo- 
ing mergers and elliptical galaxies on the basis of their 
luminosity densities, one needs to take into account po- 
tential differences in the ages of their stellar populations 
as well as any radial gradients in the average age within 
the galaxies. This is not possible in detail here, as we 
have no estimates available of the stellar population ages 
in the Toomre sequence galaxies. However, a simple es- 
timate can be made on the basis of the fact that the dy- 
namical merging age of several of the late-stage Toomre 
sequence gala xies has been estimated to be o f the order 
of lGyr (e.g.. iHibbard fe van Gorkomlll996l ). A stellar 
population fades in the Jf-band by ~ 1.5 mag be tween 
ages of 10 9 and 10 10 yrs (|Bruzual fc Charlot]|2003l ). cor- 
responding to 0.6 dex in luminosity. As the nuclei likely 
contain a mix of young and old stellar populations, this 
represents an upper limit to the amount of fading the 
nuclei may experience. Such significant fading would be 
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sufficient to explain much of the excess luminosity den- 
sity compared to E/SO galaxies for most of the Toomre 
sequence nuclei (see Figure 1251) . Also, simulations show 
that the young stars formed in mergers are likely to be 
more concentrated toward the center than the underly- 
ing older stars. A fading of the young population with 
time will therefore decrease not only the amplitude but 
also the power-law slope of the total projected luminosity 
density profile, moving the merger remnants toward the 
locus of E/SO galaxies in Figurel2"4l The effects combined 
suggest that, if left to evolve for several more Gyrs, it is 
quite possible that the properties of the Toomre sequence 
nuclei would look more similar to the nuclei of normal 
E/SO galaxies. The Toomre sequence remnants would 
then populate the steep power-law, high luminosity den- 
sity ends of the E/SO distributions shown in Figures |2"41 
andl2"5l respectively. 

One possible surprise is that the observed surface 
brightness profiles are smooth and do not show any ob- 
vious kink or upturn. This differs f rom the profiles ob- 
tained in the merge r simulations of IMihos fc HernquisH 
(|1994) and later by iSpringell (|2000D . which had a pro- 
nounced "cusp" inward of a few hundred parsecs due to 
the presence of the young starburst population (although 
see ISpringel et al.l (|2005f ) for a discussion of how black 
hole feedback may moderate these effects). Most of the 
surface brightness profiles presented here cover radii in 
the 10 pc to 1 kpc range, so we would have detected any 
such cusps if they existed. However, it should be noted 
that pronounced kinks are not necessarily a generic fea- 
ture of the merger simulations. They can be artifacts of 
the limited resolution of the simulations, and their pos- 
sible presence depends sensitively on the adopted (very 
uncertain) prescriptions for star formation and feedback. 
In any case, our observations indicate that if there is in- 
deed a centrally concentrated population of young stars, 
then either its brightness profile connects seamlessly with 
that of the underlying older stars, or it is more ex- 
ten ded (beyond the kpc rang e) than was predicted in 
the IMihos fe Hernquistl (|1994f l models. 

A predominant theory for the formation of cores in 
elliptical galaxies is through the evolution of supermas- 
sive black holes involved in a merger. If the two merg- 
ing galaxies both host a central black hole, then these 
black holes will sink to the center of the merger rem- 
nant through dynamical friction and form a binary sys- 
tem that interacts with the surrounding stars. These 
interactions lead to a hardening of the binary and an 
ejection of s tars from the galaxy center, th us excavating 
a core (e.g., [Milosavlic vic fc Merrittj|200l[ h The forma- 
tion of the core happens quickly, in 10 6 — 10 7 years, so 
that in this scenario one might expect to find cores in 
late-stage merger remnants. In fact, there is little evi- 
dence for cores in the latest-stage mergers of the Toomre 
sequence. However, our findings do not rule out this sce- 
nario, for several reasons: (a) the original merging part- 
ners may not both have had black holes; (b) a core may 
be present, but may not be of sufficient size to resolve 
with our observations; (c) the time scale for black hole bi- 
nary coalescence may have been shortened substantially 
by th e presence of copious amounts of gas (jEscala et al.1 
2005) so that there was little time for the scattering of 
stars and the formation of a core; (d) a core may have 
formed, but was subsequently filled in by stars formed 



from gas driven to the galaxy center during the merger; 
or (e) perhaps some minimum galaxy mass is essential for 
the formation of a core and the Toomre systems are not 
massive enough. Note with respect to the last point that 
galaxy luminosity may not be a good proxy for galaxy 
mass. The Toomre sequence galaxies do occupy the same 
range of absolute magnitudes Mk in which E/SO galaxies 
are sometimes found to have core profiles. Nonetheless, if 
the Toomre sequence galaxies have younger populations 
on a galaxy- wide scale (and not just on the nuclear scale 
discussed above), the fading of these young populations 
would make the galaxy fainter with time. One should 
then compare to E/SO galaxies that are fainter in Mk, 
which have core-type profile much more rarely. 

6. SUMMARY 

We have imaged the central regions of the 11 interact- 
ing and merging galaxies of the Toomre sequence in J, 
H, and K, using the NIC2 camera on fiST/NICMOS, 
and augmented these data with HST archival data in 
similar bands for those galaxies with existing NIR data. 
Compared to optical data, NIR data have the advan- 
tage that the images are much less affected by dust. 
The observations therefore offer a relatively unimpeded 
view into the nuclei, especially when compared to the 
results of our previous HST/WFPC2 study presented in 
Paper I. With our previous optical data it was difficult 
or impossible in several cases to identify the true galaxy 
nucleus /nuclei; by contrast, the near-IR data presented 
here generally reveal the nuclei clearly. In most cases 
we confirm the galaxy nucleus/nuclei reported from data 
in other wavebands, such as radio, X-ray or millimeter 
(CO) wavelengths. The nuclear positions identified from 
the different observations are generally consistent with 
each other. In NGC 7764A we detect a double nucleus for 
the first time, with a separation of 0'.'43 (260 pc), similar 
to the separation of the two nuclei in the merging galaxy 
Arp 220. The presence of a double nucleus in NGC 7764A 
is consistent with Toomre's initial assumption of the sys- 
tem being in the final merging stage prior to coalescence. 
In NGC 3256 S we find a highly extincted near-IR source 
coincident with the secondary nucleus previously identi- 
fied on the basis of radio and X-ray data. The extinction- 
corrected if-band luminosity of this source is similar to 
that of the nuclei in other Toomre sequence galaxies. 

Aperture photometry of the Toomre sequence nuclei 
reveals that their colors are consistent with the observed 
colors of elliptical galaxies and the predicted colors of in- 
termediate age or old stellar populations, but often with 
considerable extinction. There is a marginal trend for 
the nuclei to become bluer as a function of the merger 
stage in the Toomre sequence. Because stellar popula- 
tion differences alone cannot account for the range of 
colors we observe in the nuclei, we attribute this trend 
to indicate a dispersal of dust during the late stages of 
merging. This is qualitatively consistent with the pre- 
dictions of models of merging galaxies with black holes 
in which AGN feedback is invoked to disperse obscurin g 
dust as the merger progresses (e.g.. lHopkins et al1 F2005). 
It will be valuable to construct numerical models that 
make quantitative predictions for the dust and gas con- 
tent of interacting and merging galaxies as a function of 
merger (st)age, for comparison to our data. 

There is also a trend for the Toomre sequence nuclei to 



12 



Rossa et al. 



become more luminous (within an aperture of fixed phys- 
ical size, and after correction for extinction) as a function 
of the merger stage along the sequence. As long as there 
are no systematic errors in the analysis which depend on 
the merger stage, this trend is statistically significant at 
the 2a level. This result suggests that, as the merger 
proceeds, either the stellar density must increase, or the 
average stellar age must decrease (or both). As a cau- 
tionary note we remark that this result is based on the 
assumed merger stage, which may not necessarily reflect 
the true merger (st)age, given how the Toomre sequence 
was initially defined in the first place. In principle it 
is possible to test these scenarios in more detail using 
population synthesis models to fit the broad-band col- 
ors. However, the results of such an analysis are likely to 
be quite uncertain given the large amounts of extinction 
and the uncertainties in the dust geometry. We there- 
fore postpone a discussion of the stellar populations to 
a subsequent paper that will present HST/STIS spectra, 
obtained in the context of our investigation. 

We have calculated radial V— K color gradients around 
all the nuclei in our sample. The latest two merger rem- 
nants have rather shallow color gradients, whereas most 
other nuclei have steeper gradients. The nuclei with the 
strongest color gradients are mostly found in the mid- 
stage mergers. Several of the nuclei become considerably 
redder inward of about 5" from the nucleus, indicating 
strong concentrations of dust in the nuclei. 

Finally, we have derived the if-band surface bright- 
ness profiles for those Toomre sequence nuclei for which 
the morphology allows a meaningful isophotal analysis 
(12 out of the total 18). The profiles were fit with a so- 
called "Nuker law" to facilitate a comparison with other 
samples of galaxies observed with HST. The majority of 
the nuclei have steep profiles that can be characterized 
as power-law profiles. In general, the Toomre sequence 
galaxies tend to have steeper profiles and higher central 
luminosity surface densities than E/SO's. These find- 
ings can both be qualitatively explained if the Toomre 
sequence galaxies have young stellar populations con- 
centrated toward their centers. Such populations arc 
expected on the basis of TV-body simulations of spiral 
galaxy mergers, in which gas flows toward the center of 
the merger remnant trigger nuclear starbursts. This pro- 
cess may also be responsible for the observed trend in the 
Toomre sequence of increasing nuclear luminosity with 
advancing merger stage. If left to evolve and fade for 
several more Gyrs, it is quite possible that the proper- 
ties of the Toomre sequence nuclei would look similar to 
the nuclei of normal E/SO galaxies. Our results therefore 
support the view that mergers of spiral galaxies lead to 
the formation of early-type galaxies. 
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53.5" 12 h l m 53.0 B 53.5" 12 h l m 53.0 B 53.5 H 12 h l m 53.0 B 

R.A. (J2000.0) 

Fig. 1. — V-band (left), i-T-band (center), and V — K (right) images of the circumnuclear region in NGC4038. The nucleus is marked 
by a black square in the V-band image, as in this particular case several nuclear components are visible. The displayed field of view (fov) 
translates to a physical size of 1.3 kpc X 1.3 kpc. To visually indicate the scale we show an aperture of diameter 100 pc (solid circle) and, 
for all nuclei other than the Antennae, also a 1 kpc (dashed circle) at the bottom left of the V-band image. 
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12.0" 11.5" 12 h 46 m 11.0' 12.0" 11.5" 12 h 46 m 11.0 B 12.0 s 11.5* 12 h 46 m 11.0 B 

R.A. (J2000.0) 

Fig. 4.— As FigureQ] but for NGC4676 B (fov: 8.1 kpc x 8.1 kpc). 
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22.5' 23 h 18 m 22.0 B 22.5" 23 h lB m 22.0 B 22.5 B 23 h 18 m 22.0 B 

R.A. (J2000.0) 

Fig. 5.— As Figured] but for NGC 7592 (fov: 8.9kpc x 8.9kpc). 
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Fig. 7. — Enlarged view of the _R"-band image of the circumnuclear region in NGC 7764A, shown in inverted grey-scale and revealing the 
two nuclei. The displayed field of view translates to a physical size of 2.7 kpc X 2.7 kpc. 
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56.0" 55.5 s 18 h 12 m 55.0 B 54.5 s 56.0" 55.5" 18 h 12 m 55.0 s 54.5 s 56.0 s 55.5 s 18 h 12 m 55.0 B 54.5 s 

R.A. (J2000.0) 

Fig. 8.— As Figured] but for NGC6621 (fov: 3.9kpc x 3.9kpc). 
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R.A. (J2000.0) 

Fig. 9.— As FigurcQ] but for NGC 6622 (fov: 7.8kpc x 7.8kpc). 
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24.0" 23.5" ll h 4 m 23.0 B 24.0" 23.5" ll h 4 m 23.0 B 24.0" 23.5" ll h 4 m 23. 

R.A. (JZ0OO.O) 

Fig. 10.— As FigureQ] but for NGC3509 (fov: 9.5 kpc x 9.5 kpc). 
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35.5" 35.0" l h 24 m 34.5 B 35.5" 35.0 s l h 24 m 34.5 B 

R.A. (J2000.0) 

Fig. 11— As Figure[T] but for NGC 520 S (fov: 2.8kpc x 2.8kpc). 
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33.5' l h 24 m 33.0' 33.5" l h 24 m 33.0 B 33.5" l h 24 m 33.0 B 

R.A. (J2000.0) 

Fig. 12.— As Figure[T] but for NGC 520 N (fov: 2.8 kpc x 2.8kpc). 
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R.A. (J2000.0) 

Fig. 13.— As FigureQ] but for NGC2623 (fov: 6.8 kpc x 6.8 kpc). 
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R.A. (J2000.0) 

Fig. 14.— As FigureQ] but for NGC3256 (fov: 3.4 kpc x 3.4 kpc). 
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8" ll h 51 m 7 B 8" ll h 51 m 7 H 8" U h 51 m 7 B 

R.A. (J2000.0) 

Fig. 15.— As Figure[T] but for NGC3921 (fov: 7.2 kpc x 7.2 kpc). 
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45.0" 22 h 20 m 44.5" 45.0" 22 h 20 m 44.5 B 

R.A. (J2000.0) 

Fig. 16.— As FigureQ] but for NGC 7252 (fov: 5.8 kpc x 5.8 kpc). 
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Fig. 17. — J — H vs. H — K diagra m of the Toom re nuclei measured within a 100 pc aperture (indicated by red symbols). The sample of 
elliptical galaxies studied by Silva & Bothun (1998) is also shown (blue open squares). In addition, the Bruzual-Charlot cluster evolution 
track for a Chabrier IMF and solar metallicity is overplotted. The black symbols indicate clust er ages of 10 7 , 10 8 , 10 9 and 10 10 yrs. 
The indicated reddening vector was calculated from the interstellar extinction law values given bv lRieke feTeb ofskv (1983). The random 
error bars on our measurements are generally smaller than the symbol sizes. The actual errors are dominated by systematic calibration 
uncertainties (< 0.05 mag, see §[22J. Note that the symbols for the two late-stage mergers NGC 3921 and NGC 7252 (at H — K as 0.26 
and J — H = 0.74) are partially overlapping. 
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Fig. 18. — J — H vs. H — K diagram of the Toomre nuclei measured within a 100 pc aperture. The merger sample studied by 
Bushouse & Stanford (1992) is overplotted (open symbols). Note that we hav e plotted this diagram only to J — H = 1.7, omitting the 
two nuclei with the reddest colors (NGC 520 S and NGC 2623; see Figure [T7t . The indicated reddening vector was calculated from the 
interstellar extinction law values given by IRieke fc Lcbofsky (1985j). The random error bars on our measurements are g enerally smaller 
than the symbol sizes. The actual errors are dominated by systematic calibration uncertainties (< 0.05 mag, see § 12. 21 . Note that the 
symbols for the two late-stage mergers NGC 3921 and NGC 7252 (at H — K ?s 0.26 and J — H = 0.74) are partially overlapping. 
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Merger Stage 



Fig. 19. — J — K color measured within 100 pc (solid symbols) and 1 kpc (open symbols) apertures as a function of the Toomre sequence 
merger stage. The random uncertainties in the colors are generally smaller than the symb ol sizes. The true uncertainties in the colors are 
dominated by systematic uncertainties in the absolute photometric calibration (see § 12.21 1. The dotted horizontal line indicates the J — K 
color for a 10 10 yr population of solar metallicity (Bruzual & Chariot 2003), which provides a reasonable fit to the actual colors measured 
for E/SOs. Note that for NGC3256 there are no J-band measurements available. 
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Merger Stage 



Fig. 20. — V — K color measured within 100 pc (solid symbols) and 1 kpc (open symbols) apertures as a function of the Toomre sequence 
merger stage. The random uncertainties in the colors are generally smaller than the symbol sizes. The true uncertainties in the colors are 
often dominated by systematic uncertainties in the absolute photometric calibration (see § 12.21 1. The dotted horizontal line indicates the 
K color for a 10 10 yr population of solar metallicity (Bruzual & Chariot 2003), which provides a reasonable fit to the actual colors 
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measured for E/SOs. 
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Fig. 21. — A"-band luminosity (corrected for dust extinction) within a 100 pc (solid symbols) and 1 kpc (open symbols) aperture as 
a function of merger stage. The luminosities were corrected for dust extinction based on the observed J — K color, as described in the 
text. The random uncertainties in the luminosities are generally smaller than the symbol sizes. The true uncertainties are dominated by 
uncertainties in the galaxy distances. The solid lines are linear least-squares fits intended to guide the eye. There is a trend of increasing 
nuclear A"-band luminosity as a function of merger stage. 
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Fig. 22.- 
The results 
corrections 
0'.'2. The x- 



— Radial V — K color gradients for the regions around each of the Toomrc sequence nuclei, measured using circular annuli. 
for different nuclei in the same interacting system are shown using different symbols in one panel, as labeled. Note that no 
were made for differences in the PSF FWHM between the V and K bands, leading to spurious artifacts at radii smaller than 
axis measures radius both in arcscconds at the bottom and kiloparsecs at the top. 
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Fig. 22. — (continued) 



10 

If 8 

1 6 
« 4 
> 9 



Radius [kpc] 
0.5 1 



1.5 
— I — 1 



NGC 2623 











j i i i i i i i_ 



10 

1? 8 
£ 6 
SsS 4 

> 2 



0.5 
' — I — 1 



1.5 
I — ' 



NGC 3921 



as0Hx>oooo o o o 







12 3 4 
Radius [arcsec] 



38 



Rossa et al. 



N 


o 


O 




s 


J— 


B 




o 


+ _ 


o> 




o 
E 




CO . 








Sh 


Elliptici 


0.10.2O3 


Angle 


100 150 
■■■■■■■ ■■■■ 


OS. 


8 r 


£L 






Semimajor axis length (pc) 

100 10 



100 



1000 



NGC 4039 4 * * 



-H-H- 



O O O OO 




a 

£ 5 

1 3 

I 3 

n 5 

£ £ 

5 8 



NGC 4676 B 

I I I I I l| 1 1 I I ■ I ■ 1 1 1 h 



O O O o 





o o o a ooo caxwoaxattaoq 
i i i i 



0.2 05 

Semimajor axis length (arc sec) 

Fig. 23. — Major axis surface brightness profiles of the nuclear regions of selected Toomre sequence galaxies, where a fit was possible. 
The sub-panels show the surface brightness profiles overplotted with a Nuker model fit (solid line). The host galaxies of the individual 
nuclei are labeled. For comparison, we plot the profile of the NICMOS PSF in the first panel (star symbols, arbitrarily normalized). 
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Fig. 24. — Relation between the b est- fit power-law slope at Q ','15 (To. 15) and the total absolute i^-band magnitude. For comparison, 
the early-type E/SO galaxies from the Ravindranath ct al. (2001) sample are plotted as open squares. The random uncertainties in To. 15 
are small. However, the true uncertainties in To. 15 are dominated by systematic uncertainties that are difficult to quantify, having to do 
primarily with the accuracy of PSF deconvolution. Nonetheless, these uncertainties can be kept to a minimum by measuring T at a radius 
larger than the FWHM of HST's PSF. For the adopted radius of 0"15 (i.e. To. 15), the uncertainties are small enough so as to not affect the 
overall distribution of points in the plot. The random uncertainties in the absolute magnitudes Mk are generally smaller than the symbol 
sizes. The true uncertainties in Mg are dominated by uncertainties in the galaxy distances. 
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Fig. 25. — Relation between the luminosity surface density S and the total absolute K-band magnitude. The random uncertainties in 
£ are small. However, the true uncertainties in £ are dominated by systematic uncertainties that are difficult to quantify, having to do 
primarily with the accuracy of PSF deconvolution. Nonetheless, these uncertainties can be kept to a minimum by measuring £ at a radius 
larger than the FWHM of HST's PSF. For the adopted radius of 1 00 pc, the uncertainties are small enough so as to not affect the overall 
distribution of points in the plot. The symbols are as in Figure [24] which also discusses the uncertainties in Mx- 
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TABLE 1 

The Nuclei of the Toomre Sequence 



Seq. No. 


Galaxy 


Alt. name 


R.A. (J2000.0) 


Dec. (J2000.0) 


Comments on 


cz 


Distance 


O'.'l 


Class 








[hh mm ss.ss] 


[dd°mm'ss'.'ss] 


nuc. morph. 


[km s _1 ] 


[Mpc] 


[pc] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Early-stage 


1 


NGC 4038 


ARB 244 


12 01 53.06 


-18 52 03.6 


c,d,f 


1616 


13.8 


7 




1 


NGC 4039 


ARB 244 


12 01 53.61 


-18 53 10.7 


P 


1624 


13.8 


7 




2 


NGC 4676 A 


ARB 242 


12 46 10.09 


+30 43 54.6 


p+s, b? 


6613 


88.2 


13 


LINER 


2 


NGC 4676 B 


ARB 242 


12 46 11.21 


+30 43 22.2 


P 


6613 


88.2 


43 


LINER 


3 


NGC 7592 W 


VV731 


23 18 21.78 


-04 24 58.0 


p+r 


7280 


97.1 


47 


Sy2 


3 


NGC 7592 E 


VV731 






c,e,f 


7280 


97.1 


47 


1 


NGC 7764A E 


AM 2350-410 


23 53 23.75 


-40 48 25.2 


P+f 


9162 


122.2 


59 




4 


NGC 7764A W 


AM 2350-410 


23 53 23.72 


-40 48 25.4 


P 


9162 


122.2 


59 




Mid-stage 


5 


NGC 6621 


ARB 81 


18 12 55.34 


+68 21 48.8 


p + c,e 


6191 


84.4 


41 




5 


NGC 6622 


ARB 81 


18 12 59.40 


+68 21 15.0 


P 


6466 


84.4 


11 




6 


NGC 3509 


ARB 335 


11 04 23.63 


+04 49 43.9 


P 


7704 


102.7 


50 




7 


NGC 520 S 


ARB 157 


01 24 34.91 


+03 47 29.1 


P 


2281 


30.4 


15 




7 


NGC 520 N 


ARB 157 


01 24 33.29 


+03 48 02.6 


p?,c 


2281 


30.4 


15 




Late-stage 


8 


NGC 2623 


ARB 243 


08 38 24.14 


+25 45 16.1 


P 


5535 


73.8 


36 


LINER 


9 


NGC 3256 N 


VV65 


10 27 51.17 


-43 54 15.6 


P 


2738 


36.5 


18 




9 


NGC 3256 S 


VV65 


10 27 51.16 


-43 54 20.8 


c 


2738 


36.5 


18 




10 


NGC 3921 


ARB 224 


11 51 07.19 


+55 04 43.1 


p 


5838 


77.8 


38 


LINER 


11 


NGC 7252 


ARB 226 


22 20 44.77 


-24 40 42.6 


p.s 


4688 


62.5 


30 





Note. — Column (1) lists the Toomre sequence number. Cols. (2) and (3) list the galaxy name and alternate name. Cols. (4) and (5) give 
the coordinates of the nuclei as measured on our NICMOS F205W and F222M images. The Ict absolute accuracy of these coordinates is ~l''6. 
Col. (6) describes the nuclear morphology (as seen in the field of view of the NIC2 images, which typically extend to ~ 10 from the galaxy 
nucleus), with the following nomenclature: b — bar; c — complex; d — dusty; c — extended; f — filamentary; p — point-like (compact); r — ring; 
s — nuclear spiral. Col. (7) lists the heliocentric velocity, taken from the NASA Extragalactic Database (NED), except for NGC 4676, whose 
value was taken from the Lyon Extragalactic Database (LEDA). Col. (8) lists the distance, calculated using Hq — 75 km s~ 1 Mpc -1 . For the 
NGC 6621/22 system we have calculated the mean distance by averaging both distanc es calculated from t he velocities. For NGC 4038/4039 wc 
adopted the distance derived from a measurement of the tip of the red giant branch (j^avianc et al.M2004) . Col. (9) gives the physical scale in 
parscc corresponding to 0''l. Col. (10) lists the spectral classification for all galaxies with possible AGN contributions. These classifications were 
compiled from various literature sources, as described in Paper I. We note that this classification is by no means complete, since not all galaxies 
in the Toomre sequence have been studied in detail (or at all) spcctroscopically. The position of the true nucleus of NGC 7592 E is difficult to 
assess from our NIR images, as it has an extended morphology with two maxima, and other parts of the circumnuclcar region arc not covered by 
our aperture. The horizontal lines separate the groups that wc denote as early-stage, mid-stage, and late-stage mergers. 



iZST/NICMOS OBSERVATIONS OF THE TOOMRE SEQUENCE NUCLEI 



TABLE 2 

NICMOS OBSERVATIONS 



Galaxy 


datasct identifier 


Date 


Filter 


^ccn 


^exp 




(HST archive) 


[mm/dd/yyyy] 




H 


[s] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Mon 4038 


N8ER01 01 


06/1 2/2003 


F110W 


1.1264 


160 


NGC 4038 


N8FR01020 


06/1 2/2003 


F160W 


1.6037 


160 


NGC 4038 


N8FR01030 


06/12/2003 


F205W 


2.0658 


448 


Nnn 4039 

in vjvj i -t\j<ju 


N8ER0201 


05/26/2003 


F110W 


1.1264 


160 


NGC 4039 


N8FR02020 


05/26/2003 


F160W 


1.6037 


160 


NGC 4039 


N8FR 02030 


05/26/2003 


E205W 


2 0658 


448 


NGC 4676 A 


N8FR03010 


06/1 2/2003 


F110W 


1.1264 


192 


NGC 4676 A 


N8FR03020 


06/1 2/2003 


F160W 


1.6037 


160 


NGC 4676 A 

1 N V_X V_< '-t\J 1 "J iT. 


N8ER 03030 


06/1 2 /2003 


E205W 


2 0658 


448 


NGC 4676 R 


N8FR04010 


06/21 /2003 


F110W 


1.1264 


192 


NGC 4676 R 

IN VJVJ '-t\J 1 "J J-J 


N8ER 04020 


06/21 /2003 


El 60W 


1.6037 


160 


NGC 4676 R 


N8ER 04030 


06/21 /2003 


F205W 

J. ^ y J ' J v v 


2 0658 


448 


NGC 7592 E+W 


N8ER05010 


10/04/2002 


F110W 


1.1264 


160 


NGC 7592 E+W 


N8FR05020 


10/04/2002 


F160W 


1.6037 


160 


NGC 7592 E+W 


N8ER 05030 


10/04/2002 


E205W 


2 0658 


448 


NGC 7764A E+W 


N8FR06010 


04/24/2003 


F110W 


1.1264 


192 


NGC 7764A E+W 


N8FR 06020 


04/24/2003 


F160W 


1.6037 


192 


NGC 7764A E+W 

1 N VJV7 I t UT!iV -LJ^^ V V 


N8FR06030 


04/24/2003 


F205W 


2 0658 


448 


NGC 6621 


N8FR07010 


08/30/2002 


F110W 


1.1264 


224 


NGC 6621 


N8FR 07020 


08/30/2002 


F160W 


1.6037 


224 


NGC 6621 


N8FR 07030 


08/30/2002 


F205W 


2 0658 


448 


NGC 6622 


N8FR08010 


08/28/2002 


F110W 


1.1264 


224 


NGC 6622 


N8ER 08020 


08/28/2002 


F160W 


1.6037 


224 


NGC 6622 


N8ER 08030 


08/28/2002 


E205W 


2 0658 


448 


NGC 3509 


N8ER09010 


10/30/2002 


F110W 


1.1264 


160 


NGC 3509 


N8ER 09020 


1 /30 /2002 


El 60W 

-L 1UU V V 


1.6037 


160 


NGC 3509 


N8FR09030 


10/30/2002 


F205W 


2.0658 


448 


NGC 520 S 


N8ER1 001 


01 /18/2003 


F110W 


1.1264 


160 


NGC 520 S 


N8FR10020 


01/18/2003 


F160W 


1.6037 


160 


NGC 520 S 


N8FR10030 


01/18/2003 


F205W 


2.0658 


448 


NGC 520 N 


N8FR11010 


01 /18/2003 


F110W 


1.1264 


160 


NGC 520 N 


N8FR11020 


01/18/2003 


F160W 


1.6037 


160 


NGC 520 N 


N8FR11030 


01/18/2003 


F205W 


2.0658 


448 


NGC 2623 


N48H23010 


11/19/1997 


F110W 


1.1264 


88 


NGC 2623 


N48H23020 


11/19/1997 


F160W 


1.6037 


88 


NGC 2623 


N48H23030 


11/19/1997 


F222M 


2.2177 


120 


NGC 3256 


N4G601080 


11/28/1997 


F160W 


1.6037 


48 


NGC 3256 


N4G601090/91 


11/28/1997 


F222M 


2.2177 


80 


NGC 3921 


N49J09040 


01/10/1998 


F110W 


1.1264 


160 


NGC 3921 


N49J09020/60 


01/10/1998 


F160W 


1.6037 


192 


NGC 3921 


N49J09080/81/82 


01/10/1998 


F205W 


2.0658 


256 


NGC 7252 


N49J12040 


11/11/1997 


F110W 


1.1264 


128 


NGC 7252 


N49J12020/60 


11/11/1997 


F160W 


1.6037 


128 


NGC 7252 


N49J12080/81/82 


11/11/1997 


F205W 


2.0658 


256 



Note. — Col. (1) lists the galaxy identifier. Col. (2) lists the dataset identifier from 
the HST data archive. Col. (3) lists the date of observations. Col. (4) lists the used 
i/ST/NICMOS filter, and Col. (5) lists the central wavelength. Col. (6) lists the exposure 
time. Note that for NGC 3256 no F110W observations arc available in the archive. The 
K-band observations for NGC 2623 and NGC 3256 were performed with the F222M filter 
which covers the CO continuum (AA2.15— 2.29/i.m). 
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TABLE 3 

Positional offsets of the nuclei compared to other wavebands 



NIR-X-ray NIR-radio-cont. NIR-CO 



Galaxy AR.A. ADec A (a, 8) AR.A. ADec. A(o,<5) AR.A. ADec. A (a, 8) 

["] ["] ["] ["] ["] ["] ["] ["] ["] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



NGC 4038 


+1.05" 


+0.42 a 


1.8 


+0.609 


+1.309 


1.6 








NGC4039 


+0.11" 


-0.30 a 


1.8 


+1.35 9 


+0.30 9 


1.6 








NGC 4676 A 


0.00° 


-1.20 b 


1.8 


+0.24 h 


-0.67 h 


1.6 


+0.75° 


-1.36° 


1.7 


NGC 4676 B 


-0.30 6 


+0.20° 


1.8 


-1.50' 1 


-0.86 h 


1.6 


-0.45° 


-1.78° 


1.7 


NGC 7592 W 








+0.90* 


+0.73* 


1.6 


-0.75 p 


+1.20 p 


1.6 


NGC 7592 E 




















NGC 7764A E 




















NGC 7764A W 




















NGC 6621 








+0.60 J 


+0.80 J 


< 2.4 


-2.25 p 


-0.70 p 


1.6 


NGC 6622 




















NGC 3509 








+0.15 fc 


+0.80 fc 


< 1.9 








NGC 520 S 


+4.95 c 


+0.70 c 


> 35.0 


+0.15 ; 


-0.77 ; 


1.6 


+1.05° 


-0.70° 


< 1-7 


NGC 520 N 




















NGC 2623 


+0.15 d 


+2.77 d 


1.8 


+0.83 h 


-0.59 h 


1.6 


+1.05° 


-0.37 9 


< 1-7 


NGC 3256 N 


-1.20 e 


+1.40 e 


1.8 


-0.90 m 


+1.60 m 


1.6 








NGC 3256 S 


-0.90 e 


+1.40 e 


1.8 


-0.90™ 


+1.60 m 


1.6 








NGC 3921 


-0.30 f 


+1.31 f 


2.6 


+1.11' 1 


+0.09 h 


1.6 


+0.45° 


+1.57° 


< 1-7 


NGC 7252 


+0.15 f 


+1.60 / 


2.6 


+1.20" 


-2.10" 


< 2.4 


+0.75 r 


-0.44 r 


< 1.9 



Note. — Col. (1) lists the galaxy name. Cols. (2) and (3) list the positional offsets between the HST and X- 
ray detected nucleus. Col. (4) lists the uncertainty associated with this positional offset. Cols. (5) and (6) list the 
positional offsets between HST and radio continuum measurements. Col. (7) lists the uncertainty associated with this 
positional offset. Cols. (8) and (9) list the positional offsets between the HST and CO measurements. Col. (10) lists the 
uncertainty associated with this positional offset. The HST positions adopted for the nuclei arc the NICMOS positions 
listed in Table [T] except for NGC 3921 for which we used the WFPC2 position from Paper I (due to evidence for an 
~ 4" absolute astromctric error in the corresponding NICMOS datasct). The sources for the X-ray measurements 
arc Einstcin/Chandra/XMM-Ncwton, for the radio continuum measurements the VLA, and for the CO measurements 
various telescopes. The listed uncertainties associated with positional offsets are the quadrature sums of the HST 
astromctric accuracy and the astromctric accuraci es of the observatio ns in the ot her wavebands. Astrometric ac c uracie s 
were obtained from the following sources: HST, IPtak et all 020061) : Chandra, lAldcroftJ 1 )20021) ; ISchade et~aTl 1 120021) ; 
XMM-Ncwton/Einstcin, S. Snowden, priv. comm.; radio continuum and CO, the original literature sources of the 
data. See the text for a discussion of t he X-ray and radio nu cleus positions adopted for NGC 4038. Refe rences: (a) 
iZezas et al.l l| 2002bp, (MIReadl d2003ft, (c) [Fabbiano et al.| ||1992D , ( d) own measurement fr om a rchival data, (e ) Lira et al.| 
J2002D. ( fTlNolan et all J2004D . fql lNeff fc Ulvestad] 1120001), (h) iMcMahon _ et all d2002D. (i ) [Condon et all U990D, (j ) 
ICondon et all dl996l), (fc) H ibbard, unpublished. (Q lCarral et al.Nl990l). (m^lNeff et all J2003D, (n ) lHibbard et al.Hl994D , 
Col lYun &*I4 ibbard (20011), (p) Iono, priv. comm., ( a) iBryant & Scovillel J1999D , (rl lWang et"al1 ill 9921) . 



#ST/NICMOS OBSERVATIONS OF THE TOOMRE SEQUENCE NUCLEI 



TABLE 4 

NIR PHOTOMETRY OF THE TOOMRE SEQUENCE NUCLEI IN 100 PC APERTURES 



Galaxy 


J 


H 


K 


J - H 


J - K 


H - K 


V - K 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


IN Kj\^> 4Uoo 


1 ft Oft 


10. ol 


1 A IK 

14. / 


n oc; 

u.yo 


1 c; 1 
1.01 


n c;ft 
U.OD 


A ftT 

4.0 ( 


IN taU 4Uoy 


1 A ft 1 ") 
14. OZ 


1 Q ftV 


1 Q AO 
10.4Z 


U.yo 


i on 
1 .zu 


n o^ 
U.ZO 


a t^n 

4.0U 


Mpri /| ft 7ft A 
IN bt 40 10 A 


1 Q Qft 


1 ft no 
10. yz 


1 ft oc; 
10. ZD 


1 A A 
1.44 


O 1 1 
Z. 11 


n ftv 
U.D / 


ft TO 

D. / y 


IN vjL 40 ZD d 


1 1 ftfl 
1 ( .OU 


1 ft TQ 
10.(0 


1 ft KO 
10. DZ 


U.oZ 


l.Uo 


n oft 
U.ZD 


Q f^ft 

o.OO 


in^jO /oyz vv 


1 ft ftfl 
1D.DU 


1 ftf^ 
10. DO 


1 A Qft 
14. oD 


u.yo 


1 1A 
1. 1 4 


n to 
U. ( y 


A /1ft 
4.4D 


NGC 7764A 


18.25 


17.14 


16.61 


l.n 


1.64 


0.53 


6.05 


NGC6621 


17.27 


16.34 


15.86 


0.93 


1.41 


0.48 


4.27 


NGC 6622 


16.77 


15.77 


15.41 


1.00 


1.36 


0.36 


4.59 


NGC 3509 


17.90 


17.11 


16.76 


0.79 


1.14 


0.35 


3.53 


NGC 520 S 


18.18 


15.58 


14.03 


2.60 


4.15 


1.55 


7.32 


NGC 520 N 


16.99 


16.31 


15.99 


0.68 


1.00 


0.32 


2.80 


NGC 2623 


16.56 


14.68 


13.58 


1.88 


2.98 


1.10 


8.15 


NGC 3256 N 




13.08 


12.53 






0.55 


4.35 


NGC 3256 S 




16.27 


14.65 






1.62 


7.22 


NGC 3921 


14.71 


13.97 


13.70 


0.74 


1.01 


0.27 


3.39 


NGC 7252 


14.95 


14.21 


13.96 


0.74 


0.99 


0.25 


3.13 



Note. — Col. (1) lists the galaxy identifier. Cols. (2)-(4) list the J, H and K- 
band magnitudes in a 100 pc aperture. Cols. (5)-(S) lists the J — H, J—K, H — K 
and V — K colors. Note that for NGC 3256 no J-band observations were performed. 
For NGC 7592 and NGC 7764A only the primary nucleus was measured. The 
photometric random errors on the magnitudes and colors are very small for the 
ncar-IR bands (with a median of 0.01 mag) and arc not listed in the table. The 
actual errors are dominated by systematic calibration uncertainties (< 0.05 mag 
in each band, see § |2.2t . 



TABLE 5 

NIR PHOTOMETRY OF THE TOOMRE SEQUENCE NUCLEI IN 1 KPC APERTURES 



Galaxy 


J 


H 


K 


J - H 


J-K 


H - K 


V -K 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 4038 


12.53 


11.61 


10.77 


0.92 


1.76 


0.84 


3.99 


NGC 4039 


12.55 


11.77 


10.96 


0.78 


1.59 


0.81 


4.24 


NGC 4676 A 


14.85 


13.42 


12.67 


1.43 


2.18 


0.75 


5.42 


NGC 4676 B 


14.40 


13.55 


13.17 


0.85 


1.23 


0.38 


3.52 


NGC 7592 W 


14.13 


13.25 


12.51 


0.88 


1.62 


0.74 


4.21 


NGC 7764A 


14.99 


13.88 


13.22 


1.11 


1.77 


0.66 


4.55 


NGC 6621 


13.89 


12.82 


12.20 


1.07 


1.69 


0.62 


4.57 


NGC 6622 


13.75 


12.84 


12.46 


0.91 


1.29 


0.38 


3.61 


NGC 3509 


14.88 


14.02 


13.64 


0.86 


1.24 


0.38 


3.50 


NGC 520 S 


15.62 


12.56 


10.80 


3.06 


4.82 


1.76 


5.23 


NGC 520 N 


13.30 


13.18 


12.40 


0.12 


0.90 


0.78 


2.70 


NGC 2623 


13.85 


12.50 


11.66 


1.35 


2.19 


0.84 


3.99 


NGC 3256 N 




10.93 


10.16 






0.77 


3.89 


NGC 3256 S 




12.18 


11.02 






1.16 


5.05 


NGC 3921 


12.66 


11.83 


11.47 


0.83 


1.19 


0.36 


3.40 


NGC 7252 


12.54 


11.69 


11.30 


0.85 


1.24 


0.39 


3.25 



Note. — Col. (1) lists the galaxy identifier. Cols. (2)-(4) list the J, H and K- 
band magnitudes in a 1 kpc aperture. Cols. (5)-(8) lists the J — H, J — K , H — K 
and V — K colors. Note that for NGC 3256 no J-band observations were performed. 
For NGC 7592 and NGC 7764A only the primary nucleus was measured. The 
photometric random errors on the magnitudes and colors are very small (with 
a median of 0.01 mag) and are not listed in the table. The actual errors are 
dominated by systematic calibration uncertainties (< 0.05 mag in each band, see 

SB. 



Rossa et al. 



TABLE 6 

Surface Brightness Profile Fitting Parameters 



Galaxy fi b r b r b a (3 7 r .i5 Mk 

[Kmag/D"] ["] [pc] [mag] 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 



NGC4039 


13. 


.71 


0. 


547 


36. 


6 


0. 


.13 


1. 


46 


0. 


80 


1 


.10 


-23.78 


NGC4676 A 


14. 


.31 


0. 


429 


183. 


1 





.10 


1. 


01 


0. 


87 





.65 


-24.34 


NGC4676 B 


14. 


.80 


1. 


738 


743. 





4. 


.18 


1. 


86 


0. 


.65 





.94 


-24.13 


NGC 7592 W 


14. 


.07 


0. 


576 


271. 


3 


10. 


.00 


1. 


57 


1. 


.58 


1. 


.58 


-24.19 


NGC6621 


13. 


.11 


0. 


366 


146. 


6 


10 


.00 


0. 


92 


0. 


68 





.68 


-24.72 


NGC 6622 


11. 


.63 


0. 


112 


46. 


9 


2. 


.58 


1. 


31 


0. 


24 





.97 


-24.72 


NGC 3509 


14. 


.13 


0. 


355 


176. 


7 


10. 


.00 


1. 


21 


0. 


99 


0. 


.99 


-24.87 


NGC 520 N 


15. 


.04 


0. 


165 


24. 


1 


10 


.00 


0. 


40 


1. 


90 


1 


,19 


-23.73 


NGC 2623 


13. 


.10 


0. 


531 


190. 








.18 


3. 


21 


0. 


00 


1 


,12 


-24.22 


NGC 3256 N 


12 


.72 


0. 


673 


119. 


1 


1 


.79 


1. 


18 


1. 


30 


1. 


.29 


-24.72 


NGC 3921 


11. 


.16 


0. 


242 


91. 


2 


10. 


.00 


1. 


72 


1. 


12 


1 


.13 


-25.13 


NGC 7252 


13. 


.12 


0. 


517 


156. 


5 


10. 


.00 


1. 


08 


1. 


30 


1. 


.30 


-24.84 



Note. — Col. (1) lists the galaxy name. Col. (2) lists the surface brightness in the 
X-band at the break-radius (i.e., the surface brightness fit corresponding to the luminosity 
surface density £(, in equation 0). Cols. (3) and (4) list the break-radius of the surface 
brightness profile fitting in arcseconds and parsccs, respectively. Cols. (5)-(7) list the fit 
parameters. Column (8) lists the calculated power-law slope at a radius of 0''l5. Col. (9) 
lists the total absolute K-h&nd magnitudes, calculated from the apparent K-b&nd magni- 
tudes from NED. The majority of the apparent i-C-band mag nitudes are from 2MASS. The 
ones for the last four merger remnants were taken from Rothbcrg &i Joseph (2006). 



